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(54) Title: RECONSTRUCTION STABILIZER AND ACTIVE VISION 

(57) Abstract: A method for stabilizing the reconstruction of an imaged volume is presented. The method includes the steps of 
performing an analysis of the reliability of reconstruction of a radioactive-emission density distribution of the volume from radiation 
detected over a specified set of views, and defining modifications to the reconstruction process and/or data collection process to 
improve the reliability of reconstruction, in accordance with the analysis. Additionally, a method is described for performing ra- 
dioactive-emission measurements of a body structure. The method is dynamically defining further views for measurements, based on 
an analysis measurements results obtained during the current emission measurement procedure. The definition of the further views 
is performed without terminating the current measurement process. The preferred sets of views may then be applied immediately to 
the body structure. The method may be performed itcrativcly, to zoom in on a region of interest. The dynamically-defined views 
may be provided as the modifications defined by the method to improve the reliability of reconstruction. 
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RECONSTRUCTION STABILIZER AND ACTIVE VISION 



FIELD AND BACKGROUND OF THE INVENTION 

5 The present invention relates to radioactive-emission measurements of a 

volume. More particularly, the present invention relates to the accurate reconstruction 
of the volume, based on measurements from non-uniform views of a volume, and on 
the dynamic selection of views during the data acquisition process. Of particular 
interest is view selection for medical imaging and/or in conjunction with medical 

10 instruments, such as guided minimally invasive surgical instruments. 

Radionuclide imaging is one of the most important applications of radioactivity 
in medicine. Its purpose is to obtain a distribution image of a radioactively labeled 
substance, e.g., a radiopharmaceutical, within the body following administration 
thereof to a patient. Radioactive-emission imaging relies on the fact that in general, 

15 pathologies, such as malignant tumors, malfunctioning organs, and inflammations, 
display a level of activity different from that of healthy tissue. Thus, 
radiopharmaceuticals, which circulate in the blood stream, are picked up by the active 
pathologies to a different extent than by the surrounding healthy tissue; in 
consequence, the pathologies are operative as radioactive-emission sources and may 

20 be detected by radioactive-emission imaging. It will be appreciated that the pathology 
may appear as a concentrated source of high radiation, a hot region, as may be 
associated with a tumor, or as a region of low-level radiation, which is nonetheless 

■ 

above the background level, as may be associated with carcinoma. 

A reversed situation is similarly possible. Dead tissue has practically no pick 
25 up of radiopharmaceuticals, and is thus operative as a cold region. 

Thus radiopharmaceuticals may be used for identifying active pathologies as 
well as dead tissue. , 

In the discussion that follows, the term body structure is intended to include 
organs, portions of organs, a part of the body, and a whole body. The term organ 
30 target is intended to include pathological features within organs. These pathological 
features may be expressed, by radioactive-emission imaging, as any one of the 
following: 
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i. hot regions, of a radioactive emission intensity higher than the background 
level; 

ii. regions of low-level radioactive emission intensity, which is nonetheless 
above the background level; and 

iii cold regions, of a radioactive emission intensity, lower than the background 
level. 

Examples of radiopharmaceuticals include monoclonal antibodies or other 

agents, e.g., fibrinogen or fluorodeoxyglucose, tagged with a radioactive isotope, e.g., 

99M. , 67 „. 20L U ir 111- j- 123. 125. A .„ A 

technetium, gallium, thallium, indium, iodine, iodine and 

18 

fluorine, which may be administered orally or intravenously. The 
radiopharmaceuticals are designed to concentrate in the area of a tumor, and the 
uptake of such radiopharmaceuticals in the active part of a tumor, or other pathologies 
such as an inflammation, is higher and more rapid than in the tissue that neighbors the 
tumor. Thereafter, a radiation-emission-measuring-probe, which may be configured 
for extracorporeal or intracoiporeal use, is employed for locating the position of the 
active area. Another application is the detection of blood clots with 
radiopharmaceuticals such as ACUTECT from Nycomed Amersham for the detection 
of newly formed thrombosis in veins, or clots in arteries of the heart or brain, in an 
emergency or operating room. Yet other applications include radioimaging of 
myocardial infarct using agents such as radioactive anti-myosin antibodies, 
radioimaging specific cell types using radioactively tagged molecules (also known as 
molecular imaging), etc. 

The usual preferred emission for such applications is that of gamma rays, 
which emission is in the energy range of approximately 11-511 KeV. Beta radiation 
and positrons may also be detected. 

Radioactive-emission imaging is performed with a radioactive-emission- 
measuring detector, such as a room temperature, solid-state CdZnTe (CZT) detector, 
which is among the more promising that are currently available. It may be configured 
as a single-pixel or a multi-pixel detector. Alternatively, another solid-state detector 
such as CdTe, Hgl, Si, Ge, or the like, or a combination of a scintillation detector 
(such as Nal(Tl), LSO, GSO, Csl, CaF, or the like) and a photomultiplier, or another 
detector as known, may be used. 
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la - li schematically illustrate detecting units 102 and detecting blocks 101 of 
various geometries and constructions, and radioactive-emission-measuring probes 
associated with them. 

Fig. la schematically illustrates a detecting unit 102, formed as a single-pixel 
5 detector 104, for example, a room-temperature solid-state CdZnTe (CZT) detector, 
having a diameter D and a thickness t<i. Both the detector diameter D, or a diameter 
equivalent in the case of a non-circular detector, and the detector thickness x<j affect 
the detecting efficiency. The detector diameter determines the surface area on which 
radioactive emission impinges; the greater the surface area, the greater the efficiency. 
10 The detector thickness affects the stopping power of the detector. High energy 
gamma rays may go through a thin detector, and the probability of their detection 
increases with detector thickness. By itself, a single-pixel detector cannot generate an 
image; rather, all counts are distributed over the surface area of the detector. 

Fig. lb schematically illustrates the detecting unit 102 with a collimator 108, 
15 formed as a single cell of a diameter D, a length L, and a septa thickness x, attached to 
the detector 104. The collimator 108 may be, for example, of lead, tungsten or 
another material which substantially blocks gamma and beta rays. 

The collimator's geometry, and specifically, the ratio of D/L, provides the 
detecting unit 102 with a collection angle 8 analogous to a viewing angle of an optical 
20 camera. The collection angle 5 limits the radioactive-emission detection to 
substantially only that radioactive emission, which impinges on the detector 104 after 
passing through a "corridor" of the collimator 108 (although in practice, some high- 
energy gamma rays may penetrate the collimator's walls). 

Fig. lc schematically illustrates a block 101 of the detecting units 102, with 
25 the collimator 108, formed as a multi-cell collimator, of a cell diameter D. The 
collection angle 5 is defined for each of the detecting units 102 in the block, and each 
of the detecting units 102 forms a pixel in the block 101. 

Fig. Id schematically illustrates a radioactive-emission-measuring probe 100 
which comprises several detecting units 102, of different geometries and different 
30 collection angles 5, within a housing 107. 
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Figs, le - li schematically illustrate the block 101, formed as a combination 
of a scintillation detector (such as NaI(Tl), LSO, GSO, Csl, CaF, or the like), a 
collimator grid, and photomultipliers. 

As seen in Fig. le, the block 101, having proximal and distal ends 109 and 
111, respectively, vis a vis an operator (not shown), is formed of the scintillation 
detector 104, of a single pixel, and the collimators 108, to create the detecting units 
102. A plurality of photomultipliers 103 is associated with the single pixel 
scintillation detector 104, and with proper algorithms, as known, their output can 
provide a two dimensional image of the scintillations in the single pixel scintillation 
detector 104. In essence, this is an Anger camera, as known. 

The distal view 1 1 1 of the collimator grid is seen in Fig. If. 
Two optional proximal views 109 of the photomultipliers 103 are seen in Figs, 
lg and lh, as a square grid arrangement, and as an arrangement of tubes. 

An Anger camera 1 17, of the block 101 in the housing 107 is seen in Fig. li. 
In each of the cases of Figs, la - li, the geometry of the collimator 108 
determines the collection angle 5, wherein with no collimator, the collection angle 8, 
is essentially a solid angle of 4ti steradians. Thus, the collimator's geometiy affects 
both the detection efficiency and the image resolution, which are defined as follows: 

i. The detection efficiency is the ratio of measured radiation to emitted 

radiation; and 

ii. The image resolution is the capability of making distinguishable 
closely adjacent radioactive-emission organ targets, or the capability to accurately 
determine the size and shape of individual radioactive-emission organ targets. 

Naturally, it is desired to optimize both the detection efficiency and the image 
resolution. Yet, they are inversely related to each other. The detection efficiency 
increases with increasing collimator's collection angle, and the image resolution 
decreases with increasing collimator's collection angle. For example, when the ratio 
of D/L is 1/2, the collection angle 8 is substantially 2.5 steradians, so the cell views 
incident radiation within the confinement of about a 2.5-steradian sector. However, 
when the ratio of D/L is 1/12, the collection angle 8 is substantially 0.31 steradians, so 
the cell views incident radiation within the confinement of about a 0.3 1 -steradian 
sector. 
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Once the emission data is obtained, the data is processed to reconstruct the 
intensity distribution within the measured volume. The reconstruction process is 
generally complex, due to the large quantity of data which must be processed in order 
to obtain an accurate reconstruction. The following statistical model may be used to 

perform reconstruction. 

We assume an intensity distribution, /, defined over an input space £/, where U 
comprises a set of basic elements (e.g., pixels in two dimensional spaces, voxels in 
three dimensional spaces), and I(u) is the intensity of a given basic element us U. A 
detecting unit positioned on a radiation-emission-measuring-probe takes a series of 
measurements y=(vXi from different positions and orientations around the volume 
U. The geometrical and physical properties of the detecting unit, together with its 
position and orientation in a given measurement i, determine the detection probability 
400 of a photon emitted from location «. Thus the effective intensity of location u as 
viewed by the detecting unit during measurement / is <j>i(n)I(u). 

The random count A',(«) of photons that are emitted from location u and 
detected in measurement i is modeled by a Poisson process with mean <f>i(n)I(a). The 
total count of photons detected in measurement i is thus: 

yt -Poisson ( Z we u fa(u)I(u)) ( 1 a) 

or in matrix notation: 

y = Poisson(<I>I) (lb) 

where y is the vector of measurements ;>,, <& is a matrix of detection probabilities over 
measurements i and voxels i/, and I is a vector of intensity per voxel w. The 
reconstruction problem is to reconstruct the intensities I from the measurements y. 

The 2-D Radon transform is a mathematical relationship which may be used to 

reconstruct the emission intensities of volume U when the set of measurements G'Xi 
is unconstrained. The Radon transform is not statistical and does not take into account 
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the Poissonian nature of the counts. In addition, it models the views as line 
projections. The Radon transform maps the spatial domain (x,y) to the Radon domain 
(p,(|>). For a fixed projection angle, the Radon transform is simply a projection of the 
object. A technique known in the art as filtered back-projection (FBP) uses a back- 
projection operator and the inverse of the Radon transform to reconstruct the intensity 

distribution in volume U from measurements (v/)?li • 

The basic, idealized problem solved by the FBP approach is to reconstruct an 
image from its Radon transform. The Radon transform, when properly defined, has a 
well-defined inverse. However, in order to invert the transform one needs measured 
data spanning 180°. In many medical imaging situations, the positioning of the 
detecting unit relative to the emitting object is constrained, so that complete measured 
data is not available. Reconstruction based on filtered back-projection is therefore of 
limited use for medical imaging. Maximum likelihood (ML) and Maximum A 
Posteriori (MAP) estimation methods, which address the statistical nature of the 
counts, have been found to provide better image reconstructions than FBP. 

Limited-angle tomography is a reconstruction technique in the related art 
which reconstructs an image from projections acquired over a limited range of angular 
directions. The success of the reconstruction process depends upon the extent of the 
angular range acquired compared with the angular range of the missing projections. 
Any reconstruction from a limited range of projections potentially results in spatial 
distortions (artifacts) in the image. Limited angle techniques can be applied for both 
the Radon transform and the statistical models, but better results are generally 
achieved within the statistical framework. While it is known that the severity of the 
artifacts increases with the increasing angular range of the missing projections, 
limited-angle tomography does not provide information on which projections should 
be used in order to most effectively reconstruct the image. 

Maximum likelihood (ML) estimation is a widely used method in the related 
art for reconstructing an image from a constrained set of measurements. A 
parameterization of the generative model described above is obtained by assigning an 
intensity I(u) to every voxel in U, The likelihood of the observed data y = (y,),, given 
the set of parameters I = {/(*/) : w e U) is: 



Jl II • . 

• • 

• a * # 



20 
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» ^ InPobson ^«/e | ^ tfi (u) J(u)^ 

= S (-S + Vt inJ2 <M«W«) - Into!) } 



(2) 



Note that the lower and upper bound of an indexing variable (such as voxels u and 
time index t) are omitted in the following description, when they are clear from the 
5 context. 

There is currently no analytic way to solve Eqn. 2 for the maximum of the 
likelihood function. However, optimization methods that find local maxima of the 
likelihood are known. One such method is the Expectation-Maximization (EM) 
process. In EM estimation, there is no guarantee that the sequence converges to a 

10 maximum likelihood estimator. For multimodal distributions, this means that an EM 
algorithm will converge to a local maximum (or saddle point) of the observed data 
likelihood function, depending on starting values. 

Since the data generated by the model is only partially observable by our 
measurements, a basic ingredient of the Expectation-Maximization formalism is to 

1 5 define a set of random variables that completely define the data generated by the 

model. In the current case, since Yt= £u<Y<z<), the set of variables {Xu(t) : u e U; t = 
1, T} is such a set; the generated data is x = (x/) /5 where X/= (x,(u))„ ? and the 
observed data y is completely determined by x. The main tool in the EM formalism is 
the complete data likelihood: 

i*p(xin = in n i/) = e ,u n poissoii ( x ' (u > 1 / (m j) (3) 

= + .r t (ii)in(t,V(«)/(N)) + ln(.r, («)!)} 



ti 



Since the likelihood depends on the complete data, which is only partially 
observable, we take its expectation with respect to the space of the unobserved data, 
given the current set of hypothesized parameters (i.e. the current estimator). The result 
is a function Q(l\I') which assigns likelihood to sets I of model parameters, given the 

■ 

25 current set I', and given the observed data y: 



• ■ 



* • 
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Q(I|I') = E[ln P(x|I) | y; I'] (4 ) 

:-*,(i»)/(tt) +B[x c («)|»t;rjln(^ l («)/(«)) 4-C} 




t li 



where C is a term which is independent of the intensities I. The function Q(I\I r ) is 
maximized by the following new estimates: 



The expectation in Eqn. 5 is obtained as follows: 



(5) 



10 



15 



:0 



Px t {nM. ( \u)\ ltl ;i') = A . (?/t|r) 

PofasonQ/,. - xtQQI ■?'»(<')/' (<-')) PoLsson(:«,,(") I M*Q£>)) 

Poisson(7/( | Y, v <t>t(v)I(v)) 



(6) 



= Binomial (:r ( (u) j 



It follows that E[x t (u)\y ( ; I"] =>' M , and hence the EM iteration is: 



(7) 



It is provable that each EM iteration improves the likelihood. Thus, given a 
random starting estimator, the EM algorithm iterates the above improvement step 
until it converges to a local maximum of the likelihood. Several random starts 
increase the chance of finding a globally good estimator. 

It is usually desired to maximize the expected posterior probability (given a 
proper prior) rather than the expected likelihood. In that case we assume a prior 
probability on the intensities P(l) = n«*V( w ))- A proper conjugate prior for the 
Poisson distribution is the Gamma distribution: 



F(I(u)) = Gainma(I(ii) | cx u]! 3 u ) = ^ + 



(8) 



Now the maximization is done on = E[InP(x\I)p(I) |y; J^. Plugging the 

Gamma prior into g> and solving for 7(u) 9 we get the following EM iteration for the 
maximum posterior estimation: 
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/(«) = 



(9) 



(10) 



I 



The EM update step can be formulated in matrix notation as follows. Let O be 
the matrix of the projections and let I, I',y, a and J3 be represented as column 

vectors. Eqn. 10 can be written in vector and matrix notations as: 



where the explicit multiplication and division denote element-wise operations, and 
where 1 is a vector (of the appropriate length) consisting solely of l's. 

Limited computational resources (i.e., when the entire projection matrix O 
cannot be kept in memory) may require breaking the update computation according to 
a partition of <I> into a set of sub-matrices (<!>/). In that case the intensities can be 
updated gradually (using only one sub-matrix at each step) according to the following 
computation: 



where y/ is the vector of observations that are obtained using the views of <!>,-. 

In the context of image reconstruction from a constrained set of views, the 
utility of the reconstruction algorithms described above is limited. Many estimation 
algorithms, including EM and the Radon transform, require measurements from a 
complete set of views surrounding the imaged object. Although algorithms for EM 
with missing views have been developed, these algorithms are based on equally 
spaced views surrounding the imaged object, of which a number of views are not 
available. These algorithms do not provide a generalized solution for an 
unconstrained set of views, which may not be equally spaced or available for all 
directions surrounding the object. 

Singular value decomposition (SVD) is a known technique for factorizing a 
rectangular real or complex matrix, with applications in signal processing and 
statistics. SVD may be considered a generalization of Eigenvalue decomposition to 



0 + <I> r 1 




(11) 



I = 



(12) 
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m*n matrices, whereas Eigenvalue decomposition is applicable only to square 
matrices. 

SVD states that given the m-by-n matrix M whose entries are either from the 
field of real numbers or the field of complex numbers, there exists a factorization of 
the form: 



where U is an m-by-m unitary matrix, D is m-by-n with nonnegative numbers on the 
diagonal and zeros off the diagonal, and V T , the conjugate transpose of V, is an n-by- 
n unitary matrix. The elements along the diagonal of D are denoted the singular 
values. Such a factorization is called a singular- value decomposition of M. A 
common convention is to order the singular values Dj,i in non-increasing fashion, so 
that the diagonal matrix D is uniquely determined by M. 

The condition number of a matrix is defined as the ratio of the matrix's largest 
singular value to its smallest singular value. In numerical analysis, the condition 
number associated with a problem is a measure of that problem's amenability to 
digital computation. A problem with a low condition number is said to be well- 
conditioned, while a problem with a high condition number is said to be ill- 
conditioned. For example, the condition number associated with the linear equation x 
= My gives a bound on how accurate the solution y will be after approximate solution. 

SVD may be employed for the solution of linear inverse problems. The 
inverse problem for a set of linear equations, x = My, is to calculate vector y from a 
known x and M. Using SVD, with M=UDV*, the problem may be restated as x= 
UDV*y. If D is an invertible matrix, y is obtained as: 



where U and V* are unitary matrices, and D is a diagonal matrix containing the 
singular values of M. Since U and V 1 are easily transposable, a solution to Eqn. 1 3 is 
obtainable only if D is invertible. As D is a diagonal matrix of singular values ei to e n 



M = UDV T 



(13) 



y=VD" l U T x 



(14) 
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(where e k denotes the k-th singular value Dk, k ) ? O" 1 is a diagonal matrix of the 
reciprocals of the singular values. That is: 



e 



-i 



-1 



(15) 



j 



If the dimension of y is larger than the dimension of x, y cannot be determined 
directly but rather can only be estimated, which is often performed by an iterative 
procedure. If the condition number of D is very large, the multiplicative factors ef 1 
will vary greatly, thus multiplying any measurement or calculation errors based on 
individual or linear combinations of elements of y, and decreasing the likelihood of 
convergence of the iterative estimation procedure. 

Truncated SVD is a known technique for reducing sensitivity to inaccuracies 
or noise when solving a set of linear equations. In truncated SVD, the constraints 
associated with the smaller singular values are eliminated from the estimation, hi 
terms of the inverse problem of Eqn. 14, this is accomplished for the i-th singular- 
value by setting ef 1 equal to zero in matrix D" 1 . Thus the lower valued singular values 
no longer affect the estimation process. 

In order to achieve a reconstructed image which is adequate for medical 
diagnostic and treatment purposes, a reliable, high-resolution image of the tested 
object (i. e. body structure) must be obtained. Currently, reliable reconstruction 
algorithms are available only for complete data sets, which provide coverage of the 
entire volume. Such data is generally not available during medical imaging. 
Additionally, when high-resolution detecting units are used, their efficiency is 
relatively low, and the detecting units must remain at each position for a relatively 
long time in order to achieve a high probability of detection. Since during medical 
testing, measurements are generally performed at many locations as the detecting unit 
is moved relative to the observed body structure, the testing procedure generally 
requires a long time and is physically and emotionally difficult for the patient. 



• * 
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Additionally, reconstruction is based upon a large quantity of data, and is a lengthy 
and computationally complex process. 

There is thus a widely recognized need for, and it would be highly 
advantageous to have, an apparatus, system and method devoid of the above 
5 limitations, 

SUMMARY OF THE INVENTION 

According to a first aspect of the present invention there is provided a method 

10 for stabilizing the reconstruction of an imaged volume. The method includes the steps 
of performing an analysis of the reliability of reconstruction of a radioactive-emission 
density distribution of the volume from radiation detected over a specified set of 
views, and defining modifications to the reconstruction process and/or data collection 
process to improve the reliability of reconstruction, in accordance with the analysis. 

1 5 According to a second aspect of the present invention there is provided a 

reconstruction stabilizer, for improving the reliability of reconstruction of an imaged 
volume, according to a preferred embodiment of the present invention. The 
reconstruction stabilizer includes a reliability analyzer for performing an analysis of 
the reliability of reconstruction of a radioactive-emission density distribution of the 

20 volume from radiation detected over a specified set of views, and a modifier 

associated with the reliability analyzer, for defining modifications to at least one of a 
reconstruction process and a data collection process to improve the reliability of 
reconstruction, in accordance with the analysis. 

According to a third aspect of the present invention there is provided a system 

25 for generating a three-dimensional image of volume, according to a preferred 

embodiment of the present invention. The system includes a radiological imaging 
camera comprising a plurality of detectors configured for independent movement 
during data acquisition, and configured for detecting radiation emitted from the 
volume thereby to provide radiation data, a reconstructor, configured for performing 

30 an analysis of the radiation data so as to reconstruct a three-dimensional image of the 
volume, and a reconstruction stabilizer associated with the camera and the 
reconstructor. The reconstruction stabilizer includes a reliability analyzer for 
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performing an analysis of the reliability of reconstruction of a radioactive-emission 
density distribution of the volume from radiation detected over a specified set of 
views, and a modifier associated with the reliability analyzer, for defining 
modifications to at least one of a reconstruction process and a data collection process 
5 to improve the reliability of reconstruction, in accordance with the analysis, and for 
providing the modifications to at least one of the camera and the reconstructor. 

According to a fourth aspect of the present invention there is provided a 
method of radioactive-emission measurements of a body structure, according to a 
preferred embodiment of the present invention. The method includes the steps of 

10 performing radioactive-emission measurements of the body structure, at a 

predetermined set of views, analyzing the radioactive-emission measurements, and 
dynamically defining further views for measurements, based on the analyzing. 

According to a fifth aspect of the present invention there is provided a 
measurement unit for performing radioactive-emission measurements of a body 

15 structure, according to a preferred embodiment of the present invention. The 
measurement unit includes a probe for performing the radioactive-emission 
measurements of the body structure, where the probe is controllable to perform the 
measurements at a predetermined set of views, an analysis unit for analyzing the 
radioactive-emission measurements, and a view definer for dynamically defining 

20 further views for measurements, based on the analyzing. 

The present invention successfully addresses the shortcomings of the presently 
known configurations by providing methods for stabilizing the reconstruction of an 
imaged volume and of performing radioactive-emission measurements of a body 
structure. 

25 Unless otherwise defined, all technical and scientific terms used herein have 

the same meaning as commonly understood by one of ordinary skill in the art to 
which this invention belongs. Although methods and materials similar or equivalent 
to those described herein can be used in the practice or testing of the present 
invention, suitable methods and materials are described below. In case of conflict, the 

30 patent specification, including definitions, will control. In addition, the materials, 
methods, and examples are illustrative only and not intended to be limiting. 
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Implementation of the method and system of the present invention involves 
performing or completing selected tasks or steps manually, automatically, or a 
combination thereof. Moreover, according to actual instrumentation and equipment 
of preferred embodiments of the method and system of the present invention, several 
5 selected steps could be implemented by hardware or by software on any operating 
system of any firmware or a combination thereof. For example, as hardware, selected 
steps of the invention could be implemented as a chip or a circuit. As software, 
selected steps of the invention could be implemented as a plurality of software 
instructions being executed by a computer using any suitable operating system. In 
10 any case, selected steps of the method and system of the invention could be described 
as being performed by a data processor, such as a computing platform for executing a 
plurality of instructions. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 

The invention is herein described, by way of example only, with reference to 
the accompanying drawings. With specific reference now to the drawings in detail, it 
is stressed that the particulars shown are by way of example and for purposes of 
illustrative discussion of the preferred embodiments of the present invention only, and 

20 are presented in the cause of providing what is believed to be the most useful and 
readily understood description of the principles and conceptual aspects of the 
invention. In this regard, no attempt is made to show structural details of the 
invention in more detail than is necessary for a fundamental understanding of the 
invention, the description taken with the drawings making apparent to those skilled in 

25 the art how the several forms of the invention may be embodied in practice. 

In the drawings: 

Figs, la-li show detecting units and blocks of various geometries and 
constructions and radioactive-emission-measuring probes, associated with them. 

Figs, lj and lk pictorially illustrate a view and viewing parameters associated 
30 with it, in accordance with definitions of the present invention. 

Fig. 2a is a simplified flowchart of a method for stabilizing the reconstruction 
of an imaged volume, according to a preferred embodiment of the present invention. 
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Fig. 2b is a simplified flowchart of a method for analyzing a detection 
probability matrix, according to a preferred embodiment of the present invention. 

Fig. 2c is a simplified block diagram of a reconstruction stabilizer, according 
to a preferred embodiment of the present invention. 

Fig. 2d is a simplified flowchart of a method of performing radioactive- 
emission measurements of a body structure, according to a preferred embodiment of 
the present invention. 

Fig. 3 shows an object shaped as a cylinder with a front protrusion, and having 
a high-emittance portion (hotspot). 

Fig. 4a illustrates an object having two high-emission regions of interest. 

Fig. 4b illustrates the added information provided by each of views Va to V>. 

Figs. 5a and 5b are simplified flowcharts of iterative methods of performing 
radioactive-emission measurements of a body structure, according to a first and a 
second preferred embodiment of the present invention. 

Figs. 6a and b are simplified flowcharts of methods for dynamically defining 
further views, according to a first and a second preferred embodiment of the present 
invention. 

Figs. 7a and 7b present the principles of modeling, for obtaining an optimal set 
of views, in accordance with the present invention; 

Fig. 8 shows four models of a prostate emittance model. 

Figs. 9a-9f show emittance models of a given volume to illustrate view 
selection using the separability criterion. 

Figs: 9g-9i show emittance models of a given volume to illustrate view 
selection using a weighted-combination criterion. 

Fig. 10 is a simplified flowchart of an iterative method for selecting further 
views, according to a preferred embodiment of the present invention. 

Fig. 1 1 is a simplified flowchart of a single iteration of a view selection 
method, according to a preferred embodiment of the present invention. 

Fig. 12 is a simplified flowchart of a method for dynamically defining further 
views, according to a third preferred embodiment of the present invention. 
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Fig. 13 is a simplified block diagram of measurement unit for performing 
radioactive-emission measurements of a body structure, according to a preferred 
embodiment of the present invention. 

Fig. 14 is a simplified flowchart of a method for selecting an optimal initial 
predetermined set of views of a volume to be imaged, according to a preferred 
embodiment of the present invention. 

Fig. 15 shows an example of a volume with two views. 

Fig. 16 illustrates the concept of uniform coverage of a volume. 

Fig. 17 is a simplified flowchart of method for stabilizing the reconstruction of 
an imaged volume with dynamic view definition, according to a preferred 
embodiment of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present embodiments teach providing modifications of the reconstruction 
and/or imaging processes in order to obtain a reliable reconstruction of the imaged 
volume. Specifically, the methods teach analyzing the reliability of the reconstruction 
obtainable from collected emission data over a set of views to determine 
modifications which are expected to improve reliability. The present embodiments 
further teach using radioactive-emission measurements to define views for further 
radioactive-emission measurements of a body structure, to be performed during the 
current measurement process. 

With non-uniform scanning, the amount of information available for different 
voxels is not uniform, thereby constraining the ability to obtain an accurate 
reconstruction of the volume. In addition, there may be a need to focus on a region or 
regions of interest, and to control reconstruction resolution and accuracy where and 
when necessary. In such cases it is important to provide features and algorithmic 
components which handle or compensate for the lack of uniformity of information, 
and are capable of focusing on features or regions of interest, while maintaining short 
and efficient data acquisition. 

In the context of medical imaging the imaged volume corresponds to a body 
structure, which may include a whole body, portion of a body, target organ and so 
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forth. Those non-limiting embodiments presented below which are directed at 
imaging a body structure are to be understood as applying to any imaged volume. 

Before explaining at least one embodiment of the invention in detail, it is to be 
understood that the invention is not limited in its application to the details of 

5 construction and the arrangement of the components set forth in the following 
description or illustrated in the drawings. The invention is capable of other 
embodiments or of being practiced or earned out in various ways. Also, it is to be 
understood that the phraseology and terminology employed herein is for the purpose 
of description and should not be regarded as limiting. 

10 Referring now to the drawings, Figs, lj and lk pictorially illustrate a view and 

viewing parameters associated with it, in accordance with definitions of the present 
invention. 

Seen in Fig. lj is a volume U, subdivided into voxels u. The volume U is 
defined in a six-degree coordinate system (x;y;z;a>;0;cr) and has a point of origin 
15 P0(x0; yO; zO; a>0; 00; ctO). A detecting unit 102 is positioned at a location and 
orientation Pl(xl; yl; zl; col; 01; al). The detecting unit 102 has a detector 104 of a 
specific detector material of a thickness x d , and a collimator 108 of a diameter D and a 
length L, so as to define a collection angle 5. The location, orientation, and collection 
angle parameters determine a three-dimensional sector S, which is the portion of the 
20 volume U that is within the detector unit's field of view. The probability of detection 
of a given voxel outside of sector S is negligible, since a photon emitted from the 
given voxel outside of sector S will have a very low probability of reaching detector 
104 through collimator 108. 

Fig. lk schematically illustrates the emission rate of the volume U, as a 
25 function of time, given that a radioactive material of a specific half-life has been 
administered at a time TO. 

A view may thus be defined as a group of nonzero probabilities of detecting a 
radioactive emission associated with all the voxels that form sector S (Fig, lj). 
A view is sometimes referred to as a projection, and the two terms are 
30 synonymous. Furthermore, a view defined over the sector S can be naturally 
extended to be defined over the group of all voxels in the volume U, by simply 
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associating a zero probability with every voxel outside the sector S. This enables 
applying mathematical operations over the entire volume U. 

The viewing parameters, which are the factors affecting the detection of 
radioactive emissions, are as follows: 

i. location and orientation parameters: 

a location and an orientation in a six-dimensional space, Pl(xl; yl; zl; ©1; 
01; el), with respect to the origin P0(x0; yO; zO; cdO; 00; ctO) of the volume U, 
in which the detecting unit 102 is positioned; 

ii. detector-unit parameters : 

• the detector unit geometry (e.g. length L, diameter D, and/or collection 
angle 8), which together with the location and orientation parameters, 
Pl(xl; yl; zl; col; 91; al) with respect to the origin P0(x0; yO; zO; a>0; 
00; aO) define the sector S; 

• the septa thickness x, which affects the probability that a photon that 
enters the collimator will reach the detector as well as crosstalk effects 
(which occur when a photon which entered a neighboring cell 
penetrates the collimator and reaches the detector), hence, the detector 
efficiency 

• the detector material, which affects the detector efficiency; and 

• the detector thickness x<j, which affects the detector's stopping power, 
hence, its efficiency; 

Hi. attenuation parameters: 

attenuation properties of all the voxels within the sector S, as they affect the 
probabilities that radioactive emissions from a specific voxel within the sector 
S will reach the detector, wherein different voxels within the sector S may 
have different attenuation properties, since several types of tissue may be 
involved; 

iv. time parameters: 

since the radiopharmaceutical decays with a specific half-life, the time ti since 
administration, and the duration of the measurement Ati, affect the number of 
emissions that occur during the radioactive-emission measurement. 
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v. radiopharmaceutical parameters: 

The half-life ti/ 2 , of the radiopharmaceutical the types of radioactive emission, 
whether gamma or beta, and the energies of the radioactive emission affect the 
probability of detection. 

Some of these viewing parameters are fixed for a particular situation. 
Specifically, the tissue attenuation parameters are given. Additionally, the time ti 
since administration of the radiopharmaceutical is generally governed by the blood 
pool radioactivity, since it is generally necessary to wait until the blood pool 
radioactivity dies out for low-level detection to be possible. For the remaining 
viewing parameters, optimization may be carried out. 

To recapitulate the problem described above, an intensity distribution / , in 
terms of radioactive emissions per seconds, is defined over the volume U, forming our 
input space. Volume U comprises a set of basic elements u (e.g., pixels in two 
dimensional spaces, voxels in three dimensional spaces), and I(u) is the intensity in a 
given basic element u <s U. A view (also denoted a projection) q>e<& is defined by the 
set of probabilities {(p(u):ueU} 9 where (p(u) is the probability of detecting a 
radioactive emission from a voxel u, as defined by viewing parameters, such as the 
physical and geometrical properties of the detecting unit, as well as the attenuation 
parameters of the viewed volume U, and the time parameters. A measurement is 
obtained by choosing a view cpe<J>, and then sampling according to the viewing 
parameters. 

As shown in Eqn. lb, the relationship between the measurements, y, and the 
emission intensities over the body structure may be represented as y = Poisson(<M). 
The reconstruction problem we are faced with is to calculate the intensity vector I 
from the measurement vector y, given a known probability matrix <£. 

Applying S VD analysis to the probability matrix 0>, we obtain: 

<D=UDVt => y=Poisson(UDVtI) (16) 

where U and V are unitary matrices, and D is a diagonal matrix containing the 
singular values of 0>. Any destabilizing singular values in D are likewise 
destabilizing to the Poisson process. Specifically, those elements of I (or linear 
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combinations of elements) for which insufficient data is available are equivalent, 
regardless of whether or not the photon emissions are modeled as Poissonian. 

The following description is directed at a non-limiting preferred embodiment 
wherein the relationship between the measurements y and intensities I are modeled as 
5 a set of linear equations, y = <1>I, for the purposes of reconstruction. This model is 
suitable for least squares optimization when photon emissions have a Gaussian 
distribution, and thus may be reasonably applied when the measured counts (i.e. y) are 
relatively high. However the embodiments are extendable to modeling the emissions 
as Poissonian. In a preferred embodiment, when a given voxel or voxels are not 
10 within Gaussian range, the y values of several low-intensity views are combined and 
set equal to the sum of all the related linear equations, in order to bring the total count 
into Gaussian range. 

Assuming intensities with a Gaussian distribution, then by applying SVD to 
the present inverse problem, I is obtained as: 

15 

I-VD'UV (17) 

As discussed above, the condition number of D may serve as an indicator of the 
stability of the reconstruction, particularly during an iterative reconstruction process. 

20 Furthermore, known techniques such truncated SVD may be used for optimization 
and reconstruction purposes. 

The present embodiments are of a method, apparatus, and system which 
analyze the reliability of the reconstruction possible by a given imaging constellation, 
and define modifications in order to improve the quality of reconstruction. The 

25 modifications may be to the reconstruction and/or data collection aspects of the 

imaging process. As described below, the modifications may use active view selection 
during imaging and/or to guide non-uniform scanning of the imaged volume. 

Reference is now made to Fig. 2a, which is a simplified flowchart of a method 
for stabilizing the reconstruction of an imaged volume, according to a preferred 

30 embodiment of the present invention. In step 1 15, an analysis of the reliability of 
reconstruction of a radioactive-emission density distribution of the volume is 
performed. The analysis is made to determine the reliability of the data which may be 
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collected by radioactive-emission measurements taken over a specified set of views. 
In step 120, modifications to improve reliability of reconstruction are defined in 
accordance with the analysis. The modifications may involve a change or adaptation 
of the reconstruction process and/or a change or adaptation of the data collection 
5 process. Both the analysis and modification aspects of the present embodiment are 
described in detail below. 

The present embodiment tailors the reconstruction and/or data collections 
processes in order to obtain a more accurate reconstruction of the imaged volume. 
Low reliability data may cause large errors in the reconstruction process. As a result 
10 the reconstruction process may be noisy and unstable, and fail to converge properly. 
The reconstructed image may contain artifacts, unsupported frequencies and other 
errors. 

The present embodiments are applicable to all stages of the imaging and/or 
reconstruction processes. The analysis may be performed prior to data collection, in 

15 order to define a scanning procedure which provides high-reliability data. During 
data collection, the scan pattern may be adapted or views may be added (e.g. active 
vision/adaptive scanning), utilizing the new information provided by the collected 
data. After data collection, modifications may be made to the reconstruction process 
in order to counteract the effects of unreliable data and non-uniformities, and to 

20 improve reconstruction accuracy and stability. 

In a first preferred embodiment, the analysis of the reliability of reconstruction 
is based upon on analysis of the singular values of the probability matrix CD>. A non- 
limiting preferred embodiment of such an analysis is presented below. 

Reference is now made to Fig. 2b, which is a simplified flowchart of a method 

25 for analyzing a detection probability matrix, according to a preferred embodiment of 
the present invention. Instep 125, a detection probability matrix is provided. As 
described above, the detection probability matrix gives the detection probability for 
each of the views over the voxels of the volume. More explicitly, in the preferred 
embodiment each view has an associated row in the detection probability matrix. The 

30 values in a given row are derived from the probability of detection of a photon 

emerging from each of the voxels of the volume, by a detector at the associated view. 
In step 130, the singular values of the detection probability matrix are calculated. In 
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the preferred embodiment, the singular values are calculated by employing S VD on 
the detection probability matrix. Performing SVD on the detection probability matrix 
enables easy implementation of stabilization methods such as truncated SVD for 
stabilizing the reconstruction process. In step 135, those singular values which may 
have a destabilizing effect on the reconstruction process are identified as destabilizing 
singular values. 

In the preferred embodiment, a measure of the reconstruction reliability 
(denoted herein the reliability measure) is determined. The reliability measure is used 
for determining if any modification of the imaging and/or reconstruction processes is 
required, or if a stable reconstruction is possible without such modification. 

Preferably, the condition number of the probability matrix is utilized as a 
reliability measure. If the condition number is satisfactory the reconstruction may be 
considered stable, with no need for modification. Preferably, the identifying of 
destabilizing singular values is performed only if the condition number is below a 
specified magnitude. 

In an alternate preferred embodiment, the information-theoretic Fisher 
information is calculated as a reliability measure from intensity distributions provided 
of the imaged volume. An intensity distribution may be obtained either as an 
emittance model created prior to the data measurement process, or as an intermediate 
result of the reconstruction process. The Fisher information is described in detail 
below, in the context of the reliability criterion for active vision. 

In the preferred embodiment, the identification of destabilizing singular values 
is performed as follows. The ratio of the largest singular value to each of the 
remaining singular values is calculated. A singular value for which the ratio is above 
a specified threshold is considered destabilizing. One may consider that if the effects 
of the destabilizing singular values upon the reconstruction are eliminated, the 
condition number of the detection probability matrix is effectively below the 
threshold. 

Preferably, unreliable voxels or combinations of voxels associated with a 
given destabilizing singular value are identified. The intensity levels of such 
unreliable voxels or combination of voxels are unsupported in the collected emission 
data, due to the destabilizing effects of the associated singular value. Therefore the 
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accurate reconstruction of such voxels and/or combinations is unlikely, and may 
destabilize the reconstruction process. The destabilization of the reconstruction by 
unsupported voxels may lead to artifacts and unsupported frequencies in the 
reconstructed intensities. These effects are not necessarily limited to the unsupported 
regions, but may extend into other portions of the volume and create artifacts there as 
well. 

The voxels and linear combinations of voxels associated with a singular may 
be determined via SVD decomposition of the intensity distribution matrix <I>. The 
matrix V obtained by SVD decomposition of <I> indicates voxels or linear 
combinations of voxels whose intensity after reconstruction is most highly affected by 
each of the singular values in D. Each row of V indicates a weighting associated with 
a given singular value for each voxel in I. Consider the following example of a 
volume with four voxels il-i4: 



y=UDV T I-U 
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In a first example, assume that e4 is a destabilizing singular value, and that the 
fourth row of V, v 4) equals [0.001 0.01 0.97 0.0 15] . Examining v 4 it is seen that 
the coefficient for the single voxel i3 is significantly larger the coefficients of the 
other voxels. It is therefore quite apparent that the reliability of data regarding voxel 
i3 is particularly low, and it is unlikely that i3 can be accurately reconstructed. In this 
example, the preferred embodiment is to perform smoothing of i3 relative to the 
neighboring voxels. Smoothing may be a practical approach where an inaccuracy 
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occurs in spatially localized regions of the volume. An exemplary criterion to 
determine which voxels should be smoothed, is to select voxels whose corresponding 
coefficient of \\ is above a certain threshold. 

In a second example, assume that e3 is a destabilizing singular value, and that 
the third row of V, v 3 , equals [0.5 - 0.5 0.5 - 0.5] . In the present example, the 
unreliability of the data relates to a combination of voxels which may be distributed 
over the volume, rather than to an independent, spatially located voxel or group of 
voxels. In the present example, the reliability of data regarding the linear 
combination 0.5il-0.5i2+0.5i3-0.5i4, as reflected by measurements y, is low, which 
may lead to instability and inaccuracies during reconstruction. Since there are no 
dominant voxels it may be less effective to perform smoothing. Instead, the preferred 
embodiment is to constrain intensities of the combination of voxel intensities between 
iterations. 

A proposed approach for applying such a constraint is as follows. Given a 
reconstructed intensity vector V, modify the intensities prior to performing the next 

* 

i 
t 

iteration as follows: 

r= >F-v k a(v u T I) (18) 

where v k is the row of V associated with a destabilizing singular number, e k , and a is 
a proportionality factor which may be related to the magnitude of e k . The effect of the 
voxel combination on the following iteration is therefore reduced as desired, where 
the reduction may be correlated with the dominance of the associated singular value. 

Preferably, constraints are added during the reconstruction process to reduce 
or eliminate the effects of the unreliable data associated with destabilizing singular 
values. 

In the preferred embodiment, the modifications include defining constraints on 
the reconstruction process. The reconstruction may be considered unreliable if one or 
more destabilizing singular values have been found, and the modifications may be 
made in order to reduce the effect of the destabilizing singular values upon the 
reconstruction. Examples of such constraints may include: 

1) Smoothing or piece wise smoothing for all or portions of the volume 
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2) Uniting voxels 

3) Subdividing voxels 

4) Adjusting the levels of one or more voxels after a reconstruction iteration 

5) Constraining the input data 

6) Performing reconstruction with differing resolutions for different portions 
of the volume 

The above exemplary constraints are now discussed in turn. 

Smoothing a reconstructed image is performed by calculating the intensity 
level of a voxel based on the intensities of surrounding voxels, in order to control the 
magnitude and rate of the fluctuations in intensities between voxels. For example, the 
intensity of a given voxel may be corrected to better reflect its value as an average of 
the surrounding voxels. Smoothing may be performed over the entire volume, or only 
on portions of the volume deemed unreliable. 

When a non-uniform scan is performed, it is possible that not all directions 
surrounding a voxel will have equal support. It may therefore be easier to distinguish 
voxel from some neighbors than from others. In the preferred embodiment, 
smoothing is performed directionally. In directional smoothing, the smoothing is 
weighted more strongly in certain directions than in others. 

In a first example, directional smoothing is used to overcome the limitations of 
unreliable data by smoothing the voxel in the weaker directions. In a second example, 
a higher resolution is available for some portions of the volume than is available in 
other portions. If it is important to obtain a reconstruction having a uniform 
resolution, directional smoothing may be applied reduce the resolution in highly 
reliable portions of the volume to the resolution obtainable elsewhere. 

A directional smoothing policy may be based on the Fisher information 
measure, which provides directionality information. The Fisher information matrix 
indicates directionality and cross-relationships between voxels, but may be difficult to 
calculate. In the preferred embodiment, a scalar Fisher information is calculated for 
groups of views, where each group is localized in a given direction. The scalar Fisher 
information measures may then be analyzed together, in order to determine the 
relationships between the different directions. 
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Low reliability voxels may be united or combined with surrounding voxels in 
order to increase their reliability. Voxels may combined by assigning the intensity 
level of a neighboring high-reliability voxel to the low-reliability voxel. Voxel 
merging may be carried out more extensively, by repeatedly merging the lowest 
reliability voxel with one of the neighboring voxels to form an aggregate of voxels, 
until no voxel is left with a reliability below a threshold. It will be appreciated that 
the aggregate voxels have a lower resolution locally, but all other regions with good 
coverage remain of high resolution, according to their coverage. 

An opposite approach for high-reliability portions of the volume is to 
subdivide high-reliability voxels into smaller voxels in order to improve resolution 
while still achieving a required reliability. For example, if there are 1000 different 
views independently covering 1 cubic cm and almost not affected by the surrounding 
voxels, then theoretically that volume can be divided up to about 1000 voxels, if the 
views create a linearly independent set of equations with a good condition number. 

A further possible modification of the reconstruction process may be to adjust 
the effect of certain unreliable voxels or linear combinations of voxels after each of 
one or more iterations of the reconstruction process. For example, the magnitude of 
unreliable voxels or combinations of voxels associated with a destabilizing singular 
value may be reduced towards zero or towards the value of neighboring voxels by a 
specified factor, possibly in proportion to their respective weightings in the associated 
row of matrix V and/or to the associated singular value. 

An additional possible modification is to constrain the value of the input. For 
example, a realistic range of counts may be defined overall or per view based on view 
parameters, typical spatial structures and the like. Examples include adding a prior 
constraint such as a gamma distribution of intensities, a constant or piecewise- 
constant progression of intensities, a linear or piecewise-linear progression of 
intensities, smoothing or piecewise smoothing constraints, an intensity distribution 
based on the shape of or magnitude of the object being scanned (possibly determined 
during a pre-scan), or a maximal or known range of expected intensities, 

■ 

The reconstruction process may be performed in a manner that obtains varying 
resolutions over the volume. In one approach, the entire volume is defined as a single 
huge voxel, and split over the reconstruction iterations to form smaller voxels. The 
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process is performed repeatedly over all or some of the voxels of the volume, as long 
as the result of the split maintains a reconstruction reliability high enough for stable 
results. 

In a further preferred embodiment, the information-theoretic Fisher 
5 information is used as a reliability measure. Modifications to the reconstruction 
and/or data collection processes may be performed when the Fisher information is 
deemed to be outside a specified range. The calculation of the Fisher information is 
described in detail below, in the context of the reliability criterion for active vision. 
The Fisher information may be calculated from the results of a previous 
10 reconstruction iteration, or, initially, from an emittance model provided of the imaged 
volume. 

The following addresses a preferred embodiment in which the modification is 
implemented by defining views for imaging the volume. The defined views serve to 
guide the data collection process in order to obtain measurements which enable 

15 performing a stable and accurate reconstruction of the intensity distribution of the 
volume. In this manner, detecting resources may be invested effectively in order to 
improve reconstruction reliability. Scanning resources include, for example, detector 
dwell time, number of detectors, angular and translational increments, and the like - 
features that increase the amount of data collection. 

20 Such modification may be performed when a region is interesting but is still 

left with coarse resolution, so as to allocate more scanning resources, such as dwell 
time, number of detectors, angular and translational increments, to cover that region 
and to form more independent views such as to increase the reliability of the 
reconstruction of that area. Preferably, the views are defined so as to obtained a 

25 desired resolution over the region of interest. 

The defined views may yield a non-uniform scanning procedure. A scanning 
density may be specified by the angular and translational increment size, or steps, the 
smaller the increment, the higher the density. Additionally, the scanning density may 
be defined by the acquisition time at each position - the longer the time, the higher the 

30 density. 

Non-uniform scanning may be defined by specifying varying scan densities 
for imaging the volume. For example, the modification may entail adjusting a local 
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scanning density to scan a region of interest with high density, and to scan other 
regions with low density. 

The non-uniform scans may relate to non-uniform angular steps of the detector 
along a sweep, non-uniform detector translational steps, or different steps by different 
5 detectors. Some detectors may employ dense steps and others may employ sparse 
steps, for example, based on active vision as taught hereinbelow. 

The scan density may be adapted to the distance to the object of interest. 
Since resolution decreases with distance, the higher density may compensate for 
increased distance. 

10 Additionally or alternatively, the angular steps may increase in density when 

scanning a region of interest, and decrease in density, when scanning other portions of 
the volume. 

Furthermore, more than one region of interest may be scanned with dense 
steps, simultaneously. The two regions of interest may be, for example, a tissue 
1 5 region and a blood pool region. This has applicability, for example, to dynamic 
studies of blood perfusion, by providing even scanning resources both to the blood 
and to the tissue. 

Additionally, convex scans may be employed. 

Variable scans, where a same region is scanned first with a first density and 
20 then with another density, may be employed. Alternatively, a same region may be 

scanned by a first group of detectors with a first density and then by a second group of 
detectors with another density, concurrently, or at different times. Thus the same 
region is scanned with at some density by a given detector and at a different density 
by another detector. 

25 In a further preferred embodiment, view definition is performed dynamically 

during radioactive-emission measurements of the volume. The definition of further 
views for measurement during the data collection process is denoted active vision 
herein and is described in detail below. Active vision may be performed 
independently or in conjunction with the present preferred embodiment of stabilizing 

30 the reconstruction of an imaged volume. 

In the preferred embodiment, the method includes the further step of 
iteratively reconstructing the radioactive-emission density distribution of the volume, 
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preferably by EM estimation. Reconstruction reliability is preferably evaluated after 
every iteration. 

Reference is now made to Fig. 2c, which is a simplified block diagram of a 
reconstruction stabilizer, according to a preferred embodiment of the present 
5 invention. Reconstruction stabilizer 140 includes reliability analyzer 145 which 

analyses the reconstruction reliability of a radioactive-emission density distribution of 
a volume. The reconstruction reliability may be determined from measures such as 
the condition number of the probability density matrix or the Fisher information, as 
described above. Modifier 150 defines modifications to improve the reliability of the 

10 reconstruction, based on the analysis performed by reliability analyzer 145. Modifier 
150 preferably applies at least one constraint to the reconstruction process and/or data 
collection. The constraints may include applying smoothing, uniting voxels, adjusting 
the levels of voxels, defining new views or non-uniform scanning, and other 
constraints described above. The views may be defined before or during emission 

15 data collection. 

In the preferred embodiment, reconstruction stabilizer 140 is included in a 
system for generating a three-dimensional image of volume. In the preferred 
embodiment, the system further includes radiological imaging camera 155 and 
reconstructor 160. Camera 155 includes detectors capable of independent movement 

20 during data acquisition, and which are capable of detecting radiation emitted from the 
volume thereby to provide radiation data. In accordance with embodiments of the 
present invention, each block of the imaging camera construction or each detecting 
unit, where single-pixel detecting units are used, may be provided with at least one, 
and preferably, two, three, or as many as six degrees of motion such as, for example, 

25 rotational motion around the x, y, and z, axis, oscillatory motion about these axes, or 
translational motion along these axes. Camera 155 thus enables flexible view 
definition by modifier 150, increasing the likelihood that a reliable reconstruction will 
be achieved. Reconstructor 160 analyzes the radiation data provided by camera 155 
so as to reconstruct a three-dimensional image of the volume. 

30 Reconstructor stabilizer 140 guides camera 155 and reconstructor 160, by 

providing one or both with the modifications determined to improve reconstruction 
reliability. 
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In practice, it is generally necessary to reconstruct volumes at a high 
resolution (i.e. a large number of voxels) from many measurements. Implementing 
the techniques described above then requires the manipulation of extremely large 
matrices. The large size of these matrices poses difficulties both for performing the 
5 calculations and for memory management. These difficulties are particularly 
problematic when active view selection is performed. 

In the preferred embodiment, calculations are performed in a localized 
manner. That is, calculations are performed for different sections of the volume in 
turn, progressing through the volume until it is covered in its entirety. When 
10 calculations are being made of a particular section of the volume, the variables 

relating to the other sections are "frozen". Thus first a selected section of the volume 
may be stabilized, so that the reconstruction of following sections is based on data 
with improved reliability. Localized calculations, such as SVD decomposition, are 
performed on smaller sub-matrices rather than on a single large matrix, alleviating the 
15 difficulties of manipulating very large matrices. 

The division of the volume into sections may be devised in any practical way, 
for example as sequential slices or by another spatial connection. The division may 
be based on the scan pattern, so that well-supported sections may be stabilized first, 
followed by less supported section. The division may also or alternately be based on 
20 knowledge of the body being imaged, such as its shape or composition. When 

imaging a body structure, the sections may be based on the known axes or structures 
of the organ 

The following embodiments are of a method for determining further views for 
the imaging of a body structure, denoted active vision herein. Active vision addresses 

25 the problem of ensuring that the quality of data gathered during the measurement 
process is adequate to provide a high quality image. The collected data and/or the 
image reconstructed from the collected data is analyzed the while the measurement 
process is taking place. Based on the analysis, further views are defined. Since each 
view is associated with known values of the viewing parameter(s), selecting a view 

30 effectively specifies known viewing parameter values. The defined further views thus 
define a set of viewing parameter values, which are used during the current 
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measurement process in order to collect data which yields a high-quality 
reconstruction of the body structure. 

Active vision may be performed independently, as described in the 
embodiments presented below. Additionally or alternatively, active vision may be 
5 performed in conjunction with the above-described method for stabilizing 

reconstruction, by modifying data collection by dynamically providing additional 
views to improve reconstruction reliability (see Fig. 17). 

The following embodiments are not confined to a specific reconstruction 
algorithm. Further views are preferably defined based on one or more of the 
10 following: 

1) Detector photon count 

2) Geometric properties of the reconstructed body structure 

3) Information theoretic measures that quantify the quality of the data fed to 
the reconstruction algorithm 

1 5 Each of these criteria is discussed in detail below. 

■ 

Reference is now made to Fig. 2d, which is a simplified flowchart of a method 
of performing radioactive-emission measurements of a body structure, according to a 
preferred embodiment of the present invention. In step 200, radioactive-emission 
measurements of the body structure are performed at predetermined views, preferably 

20 in vivo. Preferably the measurements are performed for diagnostic purposes. These 
predetermined views are selected prior to the measurement process, based on a model 
of the body structure being imaged. In the model more and less informative viewing 
directions have been identified. The predetermined views of step 200 preferably 
include those views expected to be informative, based on an analysis of the model. 

25 Preferably the body structure is all or a portion of: a prostate, a heart, a brain, a 

breast, a uterus, an ovary, a liver, a kidney, a stomach, a colon, a small intestine, an 
oral cavity, a throat, a gland, a lymph node, the skin, another body organ, a limb, a 
bone, another part of the body, and a whole body. 

In step 210 the radioactive-emission measurements are analyzed. Preferably 

30 the analysis includes one or more of: 

1) Analyzing detector photon count(s) 

2) Analyzing detector photon count rate(s) 
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3) Identifying detector saturation 

4) Reconstructing a body structure image from emission measurements 

5) Identifying geometric properties of the reconstructed image 

6) Applying information-theoretic measures to the reconstructed image 

In step 220, further views for measurements are dynamically defined, based on 
the analysis performed in step 210. Preferably, each of the views is associated with 
viewing parameters selected from the group consisting of: detector unit location, 
detector unit orientation, collection angle, and measurement duration. Defining a 
view consists of providing a value for each of the parameters associated with the 
given view. The analysis (step 210) and/or dynamic view definition (step 220) may 
take into account external parameters including: measurement duration, time elapsed 
from the administration of the pharmaceutical to the measurement, 
radiopharmaceutical half life, radioactive emission type, and radioactive emission 
energy. 

Each of these analysis techniques, and their application to view definition, is 
now discussed in turn. While each of the analysis/view determination techniques is 
discussed as a separate embodiment, multiple techniques may be used together to 
obtain the desired image quality. 

In a first preferred embodiment, a photon count analysis ensures that the 
photon count at a given view yields an acceptable measurement error. As discussed 
above, the radiative emissions of the body structure being imaged is a Poisson 
process. In a Poisson process the Poisson noise grows inversely to the square root of 
the number of photons detected. In other words, if N photons are collected from a 
given view, the resulting signal to noise ratio (SNR) equals: 



The unprocessed detector photon count at a given view thus provides 
significant information regarding the quality of the information obtained at a given 
view. If the photon count is too low, it may be desired to continue to collect photons 
at the current location/orientation in order to obtain a satisfactory SNR. Alternatively, 
it may be determined that enough photons have already been collected, and to 
terminate the current view and move on to the next view. 
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The analysis is preferably performed by defining a global or local required 
measurement error, and comparing the square root of the obtained photon count to the 
required measurement error. Photon count analysis can be applied to the current 
and/or previous views. When a photon count of a current view is found to be too low, 
5 the duration of the current view is preferably extended in order to obtain the required 
error value. When a photon count of a past view is found to be too low, an emission 
measurement at substantially the same location and orientation but having a longer 
duration than previously is preferably performed. Alternately or additionally, the 
collection angle at the given location/orientation is preferably increased. 
0 In an additional preferred embodiment, a detector photon count is analyzed to 

identify detector saturation at a given view. Preferably, when a detector is determined 
to have saturated, a new view or views are selected to reinforce those views that have 
saturated. In an alternate preferred embodiment, further views are defined to avoid 
highly-radiating portions of the body structure. 
5 In a second preferred embodiment, a photon collection rate at a given view is 

analyzed to determine if it is within a specified range. In the preferred embodiment, 
the photon count rate is used to identify regions of high or low interest. In prostate 
imaging, for example, a region of high interest may be identified by a high photon 
rate, indicative of a tumor. In a second example, a region of high interest may be 
^0 identified in heart imaging by a low photon rate, indicative of non-functional tissues. 
After one or more areas of high and/or low interest are found, further views are 
preferably defined by selecting views to concentrate on regions of high interest and/or 
to avoid regions of low interest. It is thus possible to zoom in on a suspected 
pathology without repeating the emission measurement process. 
25 In a further preferred embodiment, the analyzing of step 210 includes 

reconstructing a radioactive-emission density distribution of the body structure. 
Reconstruction may be performed according to any applicable technique known in the 
art. The reconstruction is then used as the basis for further analysis. 

Reconstruction based on the data collected from the predetermined views 
30 provides information regarding the quality of information obtained from the preceding 
measurements, and which further views are likely to be most informative. Selecting 
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new views based on reconstruction is intended to bring us into viewing from the more 
informative views or combinations of views. 

Reference is now made to Figs. 3 and 4a-4b, which pictorially illustrate how 
different views provide differing types and quality of information. Fig. 3 shows an 
object 300 shaped as a cylinder with a front protrusion, and having a high-emittance 
portion (hotspot) 310. Four views of object 300 are shown, which can be seen to 
provide different levels of information. Front views, such as Vj, provide little 
information regarding the shape of object 300 and have relatively little attenuation 
between the detector and hotspot 310. Side views, such as V 2 , provide edge 
information regarding the object shape or profile, and when correlated with front 
views help locate hotspot 310 spatially within object 300. Top views, such as V 3 , 
provide information regarding the cylinder edge region 320. Finally, rear views, such 
as V4, are uninformative about the shape of object 300 and have high attenuation 
relative to hot region 310. 

Figs. 4a and 4b demonstrate how the proper selection of views may improve 
the quality of information obtained for the body structure, for example in 
distinguishing between two regions of interest within a given volume. 

Fig. 4a illustrates an object 400 having two high-emission regions of interest 
(ROI), 410 and 420. For clarity the views V A to V> are shown as lines in Fig. 4a, 
however in practice they will each have a finite collection angle 8. The position of 
ROIs 410 and 420 are assumed to have been estimated based on a model of object 400 
and/or a previously performed prescan. The goal of the present invention is to select 
an additional new view or views which increase the information we have regarding 
the separation of ROIs 410 and 420 within object 400. 

In simple terms, consider the object as having three regions: ROI 410 with 
intensity Ij, ROI 420 with intensity I 2 , and a low-emission region 430 between the two 
ROIs with intensity I3. The detected intensity at a given detector is proportional to 



n / 2 , where I n is the emission intensity of region n and n is the distance of region n 



from detector Vj. 

Fig. 4b illustrates the added information provided by each of the shown views, 
Va to Vp. Views V B and Vc collect emissions from all three regions, and are 
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therefore least informative. Views Vq and Ve collect emissions from only low 
emittance region 430, and therefore provide most information regarding the location 
of each ROI within the volume and the separation between ROIs 410 and 420. Views 
V A and Vf pass only through a single ROI, and therefore provide an intermediate 
5 level of information. It is a goal of the present invention to determine, while the 
emission measurements of the body structure are taking place, that views in the 
vicinity of Vd and Ve are highly informative, and to add these further views to the 
measurement process. 

A body structure reconstruction can be utilized in several ways to define 

10 further views. A first way is to identify interesting portions of the contour and 

structure of the reconstruction. For example, it is seen in Fig. 3 that top views are 
informative about edge region 320. Further top view measurements will therefore be 
informative re edge region 320, and may enable defining the edge more accurately. 
In a preferred embodiment, the reconstruction is analyzed to identify textural 

1 5 edges within the reconstruction, and view definition preferably includes selecting 
views at an angle to the textural edges. In the preferred embodiment, the angle is a 
substantially sharp angle in order to provide information regarding the edge. 

In another preferred embodiment, the reconstruction is analyzed to identify 
volumetric boundaries within the reconstruction, and view definition preferably 

20 includes selecting views at an angle to the volumetric boundaries. It is expected that 
the defined views will provide information regarding the boundary and differences in 
surrounding tissues on either side of the boundary. In the preferred embodiment, the 
angle is a substantially sharp angle. 

Another way to utilize the reconstruction to define further views is to identify 

25 suspected organ targets within the reconstruction, and to select further view(s) in 

close proximity to the suspected organ targets. A suspected organ target is typically 
detected by identifying portions of the reconstruction whose emission intensity 
distribution and spatial characteristics are typical of a suspect region. 

In a first preferred embodiment, a suspected organ target is defined as a liigh- 

30 emittance portion of the reconstruction. In prostate imaging, for example, a suspected 
organ target may be a high-em ittance portion of the structure, indicating a tumor. A 
high-emittance portion is characterized by an intensity that is greater than the 
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background intensity by a factor of at least (1+a), where a is a parameter specified by 
the user. In practice, a hotspot is usually detectable only if the radiation levels within 
the hotspot are higher than the background level by a factor of 1 .5-2. a is therefore 
typically defined between 0.5-1. However, the detectability of a hotspot rises as the 
radioactive intensity of the body rises, raising the photon count. Thus, a lower value 
of a may be used when the measured body structure is in a state of high-intensity 
emittance. For example, a body structure may be characterized by relatively high 
emittance immediately following the administration of the radiopharmaceutical, and 
be characterized by lower emittance at a later time. 

In a second preferred embodiment, a suspected organ target is defined as a 
low-emittance portion of the reconstruction. In heart imaging, for example, a 
suspected organ target is defined as a low-emittance portion of the reconstruction, 
indicating non-functional tissues. A low-emittance portion is characterized by an 
intensity that is lower than the background intensity by a factor of at least (1+P), 
where p is a parameter specified by the user. 

In the preferred embodiment the further views are used immediately for 
radioactive-emission measurements. The results of the new measurements are then 
used in another analysis to define new further views for additional measurements. 
The radioactive-emission measurements may then be said to be performed iteratively. 

Reference is now made to Fig. 5a, which is a simplified flowchart of a 
iterative method of performing radioactive-emission measurements of a body 
structure, according to a first preferred embodiment of the present invention. In step 
500, radioactive-emission measurements of the body structure are performed at 
predetermined views. In step 510, an analysis is performed of the previously 
performed emission measurements. In step 520 a decision is made whether to 
continue with further measurements. If yes, in step 530 further views are defined 
based on the analysis. Subsequent iterations continue until the decision to end the 
emission measurement process. After the first iteration, the analysis performed at a 
given stage may include consideration of all or on part of the measurements 
performed during one or more previous iterations, in addition to the new 
measurements. 
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Reference is now made to Fig. 5b ? which is a simplified flowchart of a 
iterative method of performing radioactive-emission measurements of a body 
structure, according to a second preferred embodiment of the present invention. In 
the present preferred embodiment, a reconstruction of the body structure is formed in 
step 505. The analysis step 510 is then performed utilizing data provided by the 
reconstruction(s). 

Referring again to Fig. 2, preferably, analysis step 210 includes determining 
an accuracy of the reconstruction. Accuracy is preferably determined by analyzing 
the variance of the reconstructions formed over multiple iterations. Preferably, 
further views are defined in step 220 to concentrate on the region for which higher 
accuracy is required. Regions of the reconstruction having low variance provide a 
high degree of confidence regarding the accuracy of the reconstruction in the given 
region (where a portion may include the entirety of the body structure being imaged). 
Further views may be added to the current measurements until the variance is reduced 

to a required level. 

Preferably, analysis step 210 includes determining a resolution of the 
reconstruction. Resolution is preferably determined by analyzing the full width at 
half maximum (FWHM) of peak values of the reconstruction. The FWHM is given 
by the distance between points at which the reconstructions reaches half of a peak 
value. Preferably, further views are defined in step 220 to concentrate on the region 
for which higher resolution is required. 

An additional way to define future views using the reconstruction is on an 
information-theoretic basis. A quality function expressing an information theoretic 
measure is defined. The quality function rates the information that is obtainable from 
the body structure when one or more permissible views are added to current 
measurement process. Several examples of quality functions based on information- 
theoretic measures are discussed in detail below. The quality function is used to rate 
potential further views. The measurement process may then continue at those further 
views whose addition to the previous views yields a high rating. 

Reference is now made to Fig. 6a, which is a simplified flowchart of a method 
for dynamically defining further views, according to a first preferred embodiment of 
the present invention. In step 610 a quality function is provided. The quality function 
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expresses an information-theoretic measure which rates the quality of information 
obtainable from potential further views. In step 620 a set of further views is selected 
to maximize the quality function. Preferably the selected further views fulfill certain 
constraints; for example the further views may be selected from a predefined set or 
may be located in the vicinity of a region of interest within the body structure. 

In the abovedescribed reconstruction-based analyses, the quality function is 
evaluated independently for a single reconstruction of the emission intensity of the 
body structure. However, quality functions may be defined which calculate the score 
for a given set in relation to one or more reconstructions and/or emittance models. An 
emittance model is a representation of a specific radioactive-emission intensity 
distribution within the volume U, so as to model organ targets, such as hot regions, of 
a radioactive emission intensity, higher than the background level, regions of low- 
level radioactive emission intensity, which is nonetheless above the background level, 
and cold regions, of a radioactive emission intensity, lower than the background level. 
Given an object or class of objects, emittance models may be devised to reflect 
expected or typical emission patterns for the given object. 

Developing an emittance model for a particular body structure involves 
analyzing known information about the body structure to determine expected 
emission patterns of the body structure. In order to develop a model of a particular 
body structure, for example a prostate, many factors may be considered. Physical 
aspects, such as the size and shape of the prostate and the position of the prostate 
within the torso may be considered, as well as medical knowledge regarding typical 
emissions from healthy and diseased prostates. Additional information may concern 
variations between individuals, such as age, weight percentage of body fat, and the 
like. 

For simplicity, the following discussion describes the evaluation of 
information-theoretic quality functions based on emittance models only. It is to be 
understood that at least one of the emittance models is a reconstruction of the body 
structure based on past measurements. Any remaining emittance models are provided 
externally, and may be based on general medical knowledge or on information 
gathered during a previous round of emission measurements of the body structure. 
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Reference is now made to Fig. 6b, which is a simplified flowchart of a method 
for dynamically defining further views, according to a second preferred embodiment 
of the present invention. The current method differs from the method of Fig. 6a by 
the addition of steps 605-606. In step 605 a set of one or more emittance models is 
provided (where the set includes one or more reconstructions of the body structure). 
An emittance model specifies the radiative intensity of each voxel in the body 
structure. As discussed above, some of the viewing parameters affect the radiative 
intensity of the voxels in the volume, for example the type of radiopharmaceutical and 
the time since administration of the radiopharmaceutical. Therefore, the emittance 
models provided in step 605 preferably correspond to the relevant viewing 
parameters. In step 606 a collection of possible further views of the body structure is 
provided. The collection of views includes possible further views for future 
measurements, preferably based on anatomical and other constraints. Furthermore, 
the quality function provided in step 610 may utilize multiple emission models. 

Reference is now made to Figs. 7a and 7b which present the principles of 
modeling of a body structure, in accordance with the present invention. Fig. 7a 
schematically illustrates a body section 720, having a region of interest (ROI) 730. 
The region of interest 730 is associated with a body structure 740, with a specific 
radioactive-emission-density distribution, possibly suggestive of a pathological 
feature 750, termed herein an organ target 750. Additionally, there may be certain 
physical viewing constraints, associated with the region of interest 730. Model 760, 
as illustrated in Fig. 7b, represents body section 720 as a collection of voxels, having 
a specified radioactive-emission intensity distribution. 

It will be appreciated that the model 760 of the region of interest 730 may be 
based on general medical information of the body structure 740 and common 
pathological features associated with it. Additionally, the model may be based on 
information related to a specific patient, such as age, sex, weight, and body type. 
Furthermore, a structural image, such as by ultrasound or MRI, may be used for 
providing Information about the size and location of the body structure 740 in relation 
to the body section 720, for generating the model 760. 

Reference is now made to Fig. 8, which shows four models of a prostate 
emittance model. Three views are given for each of the emittance models shown. 
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Each of the views is a maximum intensity projection (MIP) along a different one of 
the three axes. For example, in an X-Y projection the intensity of a given point is 
taken as the maximum intensity for that point along a line parallel to the z-axis. Thus 
the volume in effect becomes "transparent", and the maximum intensity regions are 

5 shown clearly. Variations in emission levels are indicated by differences in shading. 
The emittance model is shaped as an ellipsoid, the typical shape of a prostate. Since a 
diseased prostate is generally characterized by one or more hot regions, each of the 
emittance models shown has a number of high-emittance portions (two hot regions in 
Models A-C and three hot regions in Model D). Each high-emittance area is an 

10 ellipsoid, with a size of approximately 1 cubic centimeter. The intensity of the 
modeled organ targets varies randomly from two to eight times brighter than the 

surrounding volume. 

In the preferred embodiment, one or more of the emittance models contains at 
least one high-emittance portion (i.e. hot region). A prostate containing a tumor, for 
1 5 example, may be modeled as an ellipsoid volume with one or more high-emittance 
portions. 

In the preferred embodiment, one or more of the emittance models contains at 
least one low-emittance portion. A diseased heart may therefore be modeled as a 
heart-shaped volume with low-emittance portions. 

20 Note that an emittance model need not contain high- or low- emittance 

portions, but may have a uniform intensity or a slowly varying intensity. 

Following are two embodiments of quality functions, both of which include 
considerations of the emission distribution of the volume, in light of the emittance 
models provided (step 605 of Fig. 6b). The information quality given by one or more 

25 further views is evaluated in relation to multiple reconstructions and/or emittance 

models of the body structure. In a first preferred embodiment the quality function is 
based on an information theoretic measure of separability. In a second preferred 
embodiment, the quality function is based on an information-theoretic Fisher 
information measure, for ensuring reliable reconstruction of a radioactive-emission 

30 density distribution. 
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In a first preferred embodiment, the quality function implements a separability 
criterion. Separability is a measure of the extent to which the measurements that are 
obtained from each pair of models can be distinguished from one another. 

The concept of separability is illustrated in Figs. 9a-9c, each of which shows 
an emittance model of a given volume 900. Each of the emittance models contains 
two high-emittance portions (hot regions), 910 and 920, which are located in 
respectively different locations in each of the models. It can be seen that the hot 
regions in emittance models 9b and 9c are in similar portions of the volume, as 
opposed to the hot regions in model 9a. It may therefore be difficult to distinguish 
between reconstructions of emittance models 9b and 9c. The separability criterion 
ensures that the selected set includes views which provide reconstructions which 
distinguish between emittance models 9b and 9c. 

Letting / be the emittance model set, a measure for the dissimilarity of any 

two given densities in J is defined. Since most state-of-the-art estimating algorithm 

are aimed at finding ML estimators, in the current example the quality function is 
based on the likelihood function. The likelihood of an estimator of /, given a set of 
Poissonian measurements y is: 



For separability, it is desired that this measure be different for each I e L 

Since the measure is a random variable that depends on the actual measurements, all 
possible pairings of emittance models should be examined to ensure that the resulting 
distributions are separable. A quality function that captures this separability is given 
by the square of the difference between the means of the distributions normalized by 
the sum of their variances: 



AO = Z - I> (")/(") + y, in £ A («)/(«) - InCy, !) 



(19) 
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(20) 



Var£(/,) + Var£(/ 2 ) 



WO 2007/010537 



PCT/IL2006/000840 



42 

The expectations and variances in Equation 20 are taken over random 
measurements y, sampled from the true intensity I\ (note that the measure is not 
symmetric). 

Since the true (unknown) intensity can be any / <= J, a projection set <S>* that 
maximizes the worst-case separability is desired. That is: 

<D*= arg max* mm JlJ2e r SEPARABILITY*^) (21) 

Scoring for separability is based on the minimum separability obtained with a 
given set of views for all of the possible pairings of emittance models from the set of 
emittance models, thereby enabling defining a desired resolution in more than one 
direction, or in more than one portion of the volume. All of the emittance models are 
modeled on a substantially identical volume. The emittance models preferably differ 
from one another in the modeled organ targets, where the modeled organ targets are 
separated by a difference of at least the required resolution (where the displacement 
which produces the required resolution is denoted delta herein). Substantially 
identical sets of views are formed from the collection of views, and each of the 
formed sets is scored with respect to each of the pairs. One of the sets of views is 
selected, based on the minimum or average score for the plurality of pairs. 

For example, assume the set of emittance models contains the three models 
9a-9c. A separability score is calculated for a given formed set of views by applying 
Equation 20 to all three pairs 9a/9b, 9a/9c, and 9b/9c. The lowest of the three 
calculated values is taken as the separability score for the formed set. Once a 
separability score has been calculated in such manner for each of the formed sets of 
views, the view set having the highest separability is selected. 

The separability criterion may be used to ensure that a required resolution is 
obtained in all or a portion of the body. In a preferred embodiment, view set selection 
for separability is performed utilizing a set of emittance models consisting of one pair 
of emittance models having substantially identical volumes but with different 
modeled organ targets. The modeled organ targets are separated by a delta in a given 
direction so as to define a required resolution in that direction and portion of the 
volume U. Substantially identical sets of views are formed from the collection of 
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views, and scored with respect to the pair of emittance models, using a quality 
function based on the separability criterion, and one of the sets of views is selected 
based on the separability scores. The selected set is thus the set which provides the 
optimum resolution in the given direction and in the vicinity of the modeled organ 
targets. 

A similar approach may be used to ensure resolution in more than one direction 
and/or portion of the volume U. Consider for example, a pair of models of 
substantially identical volumes, as follows: The model of Fig. 9d, which 
schematically illustrates the volume U, having the modeled organ target HS, whose 
center is at a location (x;y;z) H s, and the model of Fig. 9e, which schematically 
illustrates the volume U, having the modeled organ target HS\ whose center is at a 
location (x;y;z) H s>. Fig. 9e also shows modeled organ target HS of Fig. 9d 
superimposed over the present emittance model to illustrate that HS 5 is somewhat 
displaced from HS, along the x-axis, and the displacement, is denoted as deltal in the 
present example. The displacement delta may be measured, for example, in mm. 

An optimal set of views, from the standpoint of separability, is that which will 
best distinguish between the model of Fig. 9d and that of Fig. 9e. Thus, a score, in 
terms of separability is given for the pair of models, and relates to a resolution as 
defined by the difference between the models of the pair. In the present example, the 
difference is deltal along the x-axis, around the locations of HS and HS', so the score 
given by the information theoretic measure of separability will relate specifically to a 
resolution as defined by deltal along the x-axis, around the locations of HS and HS'. 
Other portions of the volume U and other directions may have different resolutions. 

For example, consider the model of Fig. 9f, which schematically illustrates the 
volume U, having the modeled organ target HS", whose center is at a location 
(x;y;z) H s", and shows HS of Fig. 9d and HS' of Fig. 9e superimposed over the present 
emittance model. HS" is displaced from HS, along the z-axis, a displacement delta2. 
and additionally, HS" is displaced from HS 9 , along the x- and z- axes, a displacement 
delta3. 

Scores, in terms of separability, may be given to all the paring combinations, 
that is the models of Figs. 9d - 9e, relating to deltal; the models of Figs. 9d -9f, 
relating to delta2, and the models of Figs. 9e - 9f, relating to delta3. An optimal set of 
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views may be selected based on its minimum or average scores for all the pairing 
combinations; for example, the optimal set of views may be that whose average score 
for all the pairing combinations is the highest. Alternatively, a weighted average may 
be applied. 

5 In a second preferred embodiment, the quality function implements a 

reliability criterion, which is a measure of how reliably the intensity distribution of a 
given object may be reconstructed from the sampled views. Since the input to the 
reconstructed algorithm is a random sample, the output estimator is also random. A 
desired property of this output is that it be reliable in the sense that similar estimators 
10 for different projected samples (i.e. different sets of measurements) of the same input 
intensity are obtained with high confidence. 

The Fisher Information, F®{I) is a measure, known in the art, which is used to 
evaluate the expected curvature of the likelihood of a model / (taken over 
measurements sampled from the model /). The Fisher Information is defined as: 

15 

F 0 (I) = -EV 2 L(I) (22) 

The derivatives are taken with respect to the parameters of the intensity /, and 
the expectation is taken with respect to the random counts. Intuitively, a sharper 

20 curvature means that a maximum-likelihood estimation algorithm is more likely to 
produce a low- variance estimator. Indeed, this property is captured by the Cramer- 
Rao lower bound, which states that the inverse of the Fisher Information is a lower 
bound for the variance of any unbiased estimator, where an estimator/of the intensity 
/ is unbiased if E/(y) = /. 

25 The Fisher information provides a value, Fq(I) UiV , for each pair of voxels u and 

v. To provide a single measure for the reliability of the estimator, the average level, 
worst case, or other reasonable measure may be taken over the voxels. In the current 
example, the quality function is based on the average Fisher information. Starting 
with a set of emittance models to be used by the estimating algorithm to provide a 

30 reliable estimator, the average of the above measure is maximized over the entire set, 
defining the selected set <D* as: 
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<D-=argmax,X IK(/)] (23) 

Alternately, a set may be chosen to minimize the inverse of the Fisher 
information, by selecting for: 

<&-=argmiti,£ ZK^'J < 24 > 

/eJT msC/ 

Since inverting F®(I) may be computationally expensive, the [/^(/y 1 J term 

in Equation 24 may be replaced with l/[Fo(/)]«.«(thus neglecting the off-diagonal 
elements of the Fisher Information matrix Note that Equations 23 and 24 are 

not mathematically equivalent, and may therefore yield different selected sets. 

In the preferred embodiment, evaluating the quality function is performed 
using the reliability criterion, and two or more emittance models are provided, having 
substantially identical volumes, but different modeled organ targets. Substantially 
identical sets of views are formed for all the emittance models, and each set is scored 
for reliability. One of the sets of views is then selected based on the average score for 
all the of the emittance models. 

In a further preferred embodiment, a weighted combination of several 
information theoretic measures is used. For example, a plurality of models may be 
provided, all having substantially identical dimensions and volumes, as follows: 

i. a pair of models with slightly different distributions of radioactive emission 
sources, as seen in Figs. 9g and 9h, for scoring sets of views in terms of separability; 

ii. a model with a given distribution of radioactive emission sources, as seen 
in any one of Figs. 9g, 9h or 9i, for scoring sets of views in tenns of reliability. 

Identical sets of views may be applied to all the models, and each view may be 
scored in terms of separability and reliability. An optimal set of views may be 
selected based on a summation of the two scores, or based on a weighted average of 
the two scores, where the relative weight given to each criterion reflects the relative 
importance of each measure per the given application. 

The quality function is preferably defined in accordance with one of the 
following: worst case effectiveness for the given further view over the volume, 
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average effectiveness for the given further view over the volume, worst case 
effectiveness for the given further view over the set of emittance models and average 
effectiveness for the given further view over the set of emittance models. 

Maximization of the quality function may be performed utilizing any method 
known in the art such as simulated annealing and gradient ascent. In the simulated 
annealing (SA) method, each point of the search space is compared to a state of some 
physical system. The quality function to be maximized is interpreted as the internal 
energy of the system in that state. Therefore the goal is to bring the system from an 
arbitrary initial state to a state with the minimum possible energy. 

The neighbors of each state and the probabilities of making a transition from 
each step to its neighboring states are specified. At each step, the SA heuristic 
probabilistically decides between moving the system to a neighboring state s' or 
staying put in state s. The probabilities are chosen so that the system ultimately tends 
to move to states of lower energy. Typically this step is repeated until the system 
reaches and acceptable energy level. 

Gradient ascent, on the other hand, is based on the observation that if a real- 
valued function F(x) 9 such as the quality function of the present embodiments, is 
defined and differentiable in a neighborhood of a point a y then F(x) increases fastest if 
one goes from a in the direction of the gradient of F at a, YF(a) . It follows that if: 

b=a+yVF(a) (25) 
for y > 0 a small enough number, then F(a)< F(b). Gradient ascent starts with a guess 
x 0 for a local maximum of F, and considers the sequence x 0 , x l# x 2 ,. . . such that: 

AW=** + >'VF(x n ),n>0, (26) 

Since F(x 0 ) < F(xj) < F(xj) < . . . , the sequence (x n ) is expected converges to a 
local maximum. 

Preferably, the set of views selected with the quality function is increased by 
at least one randomly selected view. The randomly selected view(s) increase the 
probability that the quality of information obtained with the further views is 
maximized globally rather than locally. 

As discussed above, selecting the best set of size N from amongst a large set of 
candidate projections is computationally complex. Since the size of the collection of 
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views and of the required set may be large, a brute force scheme might not be 
computationally feasible. 

In an additional preferred embodiment, a so-called "greedy algorithm" is used 
to incrementally construct larger and larger sets, until the required number of further 
5 views is defined. When multiple further views are required, it is computationally 
complex to maximize the quality function over all possible combinations of further 
views. The greedy algorithm reduces the computational burden by selecting the 
further views one at a time. The algorithm starts with a current set of views, and for 
each iteration determines a single view that yields the maximum improvement of the 
10 set score (hence the name "greedy"). 

In theoretical terms, assume p(-) is the quality measure we are using for the 
view selection, and assume without loss of generality that we are trying to maximize 
this measure. We gradually build a set W of projections as follows. We stall with an 
empty set W =0, and at every stage choose the projection that maximizes the quality 
15 measure when added to the current set: 

W<r- arg max,r (p<7'P) \W = Wu <j>€<P} (27) 

In other words, during a given iteration, a respective score is calculated for a 
20 combination of the previous set with each of the views which is not a member of the 
current set. The current set is then expanded by adding the view which yielded the 
highest respective score, and the expanded current set serves as the input to the 
following iteration. Thus the number of times the scoring function is calculated per 
iteration drops from iteration to iteration. For a large collection of possible views, the 
25 greedy algorithm reduces the total number of computations required for set selection. 

Reference is now made to Fig. 10, which is a simplified flowchart of an 
iterative "greedy" method for defining further views, according to a preferred 
embodiment of the present invention. In step 1000 a collection of views of the body 
structure is provided. The collection of views includes possible further views for 
30 future measurements, preferably based on anatomical and other constraints. In step 

1010, the set of views used for the previous emission measurements is established as a 
current set of views. In step 1020 the view set is incrementally increased by a single 



WO 2007/010537 



PCT/IL2006/000840 



48 

further view during each iteration, until the required number of further views has been 
selected. 

Reference is now made to Fig. 11, which is a simplified flowchart of a single 
iteration of the view selection method of Fig. 10, according to a preferred 
embodiment of the present invention. The method of Fig. 1 1 expands the current set 
of views by a single view. The method begins with a current set of views, which is 
the predetermined set (step 1010 above) for the first iteration of the greedy algorithm, 
or the set formed at the end of the previous iteration (step 1 120 below) for all 
subsequent iterations. In step 1 100, a respective expanded set is formed for each view 
not yet in the current set of views. A given view's expanded set contains all the views 
of the current set of views as well as the given view. In step 1 1 10, a respective score 
is calculated for each of the expanded sets using the quality function. In step 1 120, 
the view which yielded the highest-scoring expanded set is selected as a further view, 
to be used for further radioactive emission measurements. Finally, in step 1 130, the 
current set is equated to the highest-scoring expanded set by adding the selected view 
to the current set. The newly formed current set serves as an input to the subsequent 
iteration, until the desired number of views is attained. 

Reference is now made to Fig. 12, which is a simplified flowchart of a method 
for dynamically defining further views,, according to a third preferred embodiment of 
the present invention. In step 1210, a collection of possible further views for 
performing radioactive-emission measurements of the body structure is provided. 
Each of the views is associated with at least one viewing parameter. Preferably the 
viewing parameters consist of at least one the following: detector unit location, 
detector unit orientation, collection angle, and measurement duration. 

In step 1220 at least one quality function is provided. Each quality function is 
for evaluating sets of views, essentially as described above. A single quality function 
may be used to select several sets of views, where each set of views contains a 
different number of views. 

In step 1230, multiple sets of further views (where a set may include a single 
further view) are formed from the collection of views, using the quality function(s) 
provided in step 1220. In a first preferred embodiment, each of the sets is formed 
using a different one of the quality functions. In an alternate preferred embodiment, 
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one or more of the quality functions are used to form more than one set of views, 
where sets formed with the same quality function have differing numbers of views. 
In step 1240, a selected set of views is obtained from the sets formed in step 

1230. 

In a first preferred embodiment, the final set of views is obtained by choosing 
one of the sets formed in step 1230 using a set selection criterion. For example, a 
respective set is formed in step 1230 for the separability and reliability criteria 
independently. A set selection criterion which calculates an overall performance 
rating for a given set taking both criteria into account is defined, and the formed set 
with the highest overall rating is selected as the final set. 

In another preferred embodiment, the selected set of views is obtained by 
merging the sets formed in step 1230 according to the relative importance of the 
respective quality function used to form each set. 

In the preferred embodiment, the method further consists of providing at least 
one emittance model and/or reconstruction representing the radioactive-emission 
density distribution of the volume, and of evaluating with at least one of the quality 
functions of step 1220 is performed in relation to the emittance models. 

As discussed above, since each view is associated with one or more 
parameters, the selected set yields a group of parameter values for performing 
effective detection of the intensity distribution of the body structure. For example, if 
each view is associated with a view location parameter the selected set defines a set of 
locations for collecting emission data from an object, in order to provide a high- 
quality reconstruction of the intensity distribution of the body structure. 

Reference is now made to Fig. 13, which is a simplified block diagram of 
measurement unit for performing radioactive-emission measurements of a body 
structure, according to a preferred embodiment of the present invention. 
Measurement unit 1300 includes probe 1310, analyzer 1320 and view definer 1330. 
Probe 1310 performs the radioactive-emission measurements of the body structure. 
Radioactive-emission-measuring probe 1310 preferably comprises several detecting 
units, which may be of different geometries and different collection angles 8, within a 
housing (see for example Fig. Id). Preferably, the orientation and/or collection angle 
of the individual collimators is controllable. Analyzer 1320 analyzes the radioactive- 
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emission measurements obtained from probe 1310. View definer 1330 dynamically 
defines further views for measurements, based on the analysis provided by analysis 
unit 1 320. The analysis and view definition are performed substantially as described 
above. 

The abovedescribed methods for radioactive-emission measurements of a 
body structure begin by performing measurements at a predetermined set of views. 
The results of the initial measurements are then analyzed and further views are 
defined. 

The initial set of views is preferably selected based on information theoretic 
measures that quantify the quality of the data fed to the reconstruction algorithm, in 
order to obtain the best data for reconstructing a three-dimensional image of the body 
structure. 

In accordance with the present invention, our approach is delineated by the 
following process: 

i. modeling the body structure or region of interest as a model of a volume U, 
possibly with one or several modeled organ targets HS, within anatomical 
constraints AC (see Figs. 7a-7b); 

iL obtaining an optimal and permissible initial set of views for the modeled 
volume U; and 

iii. initiating the radioactive-emission measurements at the optimal set of views. 

The following section concentrates on the second step of the process, namely, 
obtaining the optimal and permissible set of initial views for performing the 
radioactive-emission measurements of the body structure. The initial predetermined 
set of views is denoted herein the optimal set of views. 

We consider here the following problem: Assume that there is a large pool of 
candidate views to choose from, but due to time restrictions or other restrictions we 
are limited to a specific number of views N. Which are the best N projections in terms 
of the quality of the reconstruction? It is further assumed that the pool of projections 
may be constrained, and hence general sampling theorems (e.g., Radon Transform) 
cannot be applied. For instance, we consider a scenario in emission tomography 
where the detecting unit can be located on top of one face of a given volume but not 
on the others. In such cases, methods such as EM estimation do not clearly establish 
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what is the best scanning scheme, in order to provide the best reconstruction of the 
radioactive intensity distribution of the volume. 

Reference is now made to Fig. 14, which is a simplified flowchart of a method 
for selecting an optimal initial predetermined set of views of a volume to be imaged, 
5 according to a preferred embodiment of the present invention. In step 1400, an 

emittance model of the body structure being imaged is provided. The body structure 
is modeled as a volume, and is preferably represented as a collection of voxels. In 
contrast to the use of emittance models for evaluation of the quality function, in the 
following all of the emittance models are provided prior to performing the 

10 radioactive-emission measurements, and thus do not include reconstructions of the 
body structure based on a current round of measurements. In step 1410, a collection 
of views of the body structure is provided. Each view reflects a possible data 
collection option, from which the goal is to select a set of views which together 
provide the desired image reconstruction quality. Each view is associated with one or 

15 more viewing parameters, such as the location and orientation of the detecting unit 
with respect to the volume, as described in detail above. 

The view set is preferably tailored to correspond to the actual limitations and 
constraints likely to be encountered when collecting data for a given object, that is a 
given body structure. For example, in medical applications the detecting unit can 

20 generally view the body structure only from certain distances and orientations, which 
are determined by the patient's anatomy and the structure of the probe. The view set 
may therefore be designed to contain views having only those viewing parameter 
values consistent with attainable distances and orientations. The view set may also be 
designed to contain views having only those viewing parameter values suitable for a 

25 given measurement scenario, for example having identical values for the type of 
radiopharmaceutical and time since administration. 

Preferably, the collection of views represents a quantized continuum of views. 
The view collection thus reflects the detection probability distribution for a detecting 
unit making periodic measurements while moving along a trajectory. 

30 Reference is now made to Fig. 1 5, which shows a non-limiting example of a 

volume 1500 with two views, Vi and V2. Volume 1500 divided into twenty-seven 
cube-shaped voxels (labeled 1-27). Each view reflects a number of viewing 
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parameters, in the case shown a location (in the XYZ coordinate system), orientation 
(co, 0 and or) and collection angle 8. The viewing parameters determine which voxels 
in volume 1500 are within a detecting unit's collection angle, and generally affect the 
probability of detection for each voxel. The probability of detection may be affected 

5 by additional factors, including the attenuation coefficients within the volume. 

During data collection, the probability of detection for each voxel is dependent 
on the parameters outlined above. In the preferred embodiment, a respective 
projection is calculated for each view, giving the view's detection probability 
distribution (i.e. the detection probability for each voxel of the volume). For a given 

10 view, the associated projection will have significant detection probabilities only for 
those voxels within the sector defined by the detecting unit's collection angle and 
location and orientation, as illustrated in conjunction with Fig. lj. 

The detection probability distribution for each view, that is the group of 
probabilities for each voxel, for a given view, is calculated according to techniques 

1 5 known in the ait for determining the detection probability of a radioactive emission 
from a given voxel, under the constraints specified by the viewing parameters, for 
example, based on computer simulations of the geometry and the other viewing 
parameters, delineated hereinabove. 

Generally, more distant voxels will have a lower probability of detection than 

20 closer voxels along the same line of sight. The volume attenuation coefficient may be 
constant or may vary over the volume. Thus, different sets of views may be produced 
for a single sector, by defining volumes with differing attenuation coefficients. For 
example, bone tissue and muscle tissue have different attenuations. Anatomical 
knowledge of the body structure being imaged may be used to develop a model of the 

25 volume U with a non-uniform attenuation that reflects the expected attenuation of the 
given body structure. 

Referring again to Fig. 14, in step 1420 a scoring function is provided. The 
scoring function rates die information that is obtainable from the volume using any set 
of views (containing at least one view) taken from the collection of views. Evaluation 

30 of the scoring function does not require data from previously performed 

measurements. This is in contrast to the quality function, which rates the information 
quality on the basis of previous emission measurements as well as by the 
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characteristics of possible further views. Preferably, the scoring function is a function 
of the detection probability distribution, as given by the projections. Several preferred 
embodiments of the scoring function, are discussed in detail below. 

In step 1425, sets of views are formed from the collection of views. A score is 
5 then calculated for each of the sets. 

In step 1430, the scores calculated in step 1425 are used to select one of the 
formed sets of views as the optimal set of views for the abovedescribed methods of 
performing radioactive-emission measurements of a body structure. A given scoring 
function may be used to select a set in a number of ways. In a first preferred 
10 embodiment, a required number of views is specified, and the highest scoring set with 
the specified number of views is selected. In a second preferred embodiment, the user 
may specify a minimal score which is known to provide satisfactory information 
quality, and select the smallest set which provides the specified score. However given 
a large collection of views the required number of calculations may be prohibitive. A 
1 5 third preferred embodiment used to reduce the computational burden is the greedy 
algorithm embodiment, similar to the abovedescribed method of Figs. 10-11. 

As discussed above, the scoring function is a measure of the quality of 
information which may be gathered for the volume using the given set of views. Both 
the separability and reliability criteria discussed above may be utilized as scoring 
20 functions, substantially as described above. 

In an additional preferred embodiment, the scoring function implements a 
uniformity criterion, to ensure uniform coverage of the volume. It is often desired to 
obtain a uniform reconstruction quality among a given set of voxels WcXJ 9 where W is 
the set of voxels for which it is desired to obtain uniform detection. Note that by 
25 selecting the appropriate W, the uniformity criterion is applied to all or a portion of 
the body. The uniformity criterion ensures that the spread of the total influence of 
each element on the set of measurements is as uniform as possible. The uniformity 
criterion depends only on the collection of views <& and requires no assumptions on 
the distribution /. 
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For a set <I> of projections, the total influence of an element, u, is given by 
£<pe<i> <P («)• Normalizing these values to P<j>(u), such that ^ P 0 (*/) = 1 , a probability 

measure is obtained for which the entropy //(O) can serve as a uniformity measure: 

The selected set O* is the set (containing the required number of views) that 
satisfies: 

<D* = arg max*ff(<X>) (29) 

Reference is now made to Fig. 16, which illustrates the concept of uniform 
coverage of a volume. Volume 1600 consists of twenty- seven cube-shaped voxels, 
including voxels 1610 and 1620. A, B, C and D are three views of volume 1600, 
showing the detector position and orientation for each view. For the uniformity 
criterion, it is desired that the overall influence of voxels 1610 and 1620 be 
approximately equal. 

Assume that the probabilities of detection are as follows: 



View 


Voxel 1610 


Voxel 1620 


A 


0.6 


0 


B 


0.2 


0.5 


C 


0 


0.3 


D 


0 


0.1 



Consider two possible sets of views: set {A, B, C} and set {B, C, D}. For set 
{A ? B, C}, the total contribution of voxel 1610 is 0.8 (0.6+0.2+0) and of voxel 1620 
is 0.8 (0+0.5+0.3). Normalizing these values for set {A, B, C} gives a probability set 
of [0.5,0.5]. For set {B, C, D} 3 the total contribution of voxel 1610 is 0.2 (0.2+0+0) 
and of voxel 1620 is 0.9 (0.5+0.3+0.1). Normalizing these values for set for set {B, 
C, D} gives a probability set of [0.18,0.82]. Thus: 
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H({A, B, C})= - (0.5 * log 2 0.5 + 0.5 * log 2 0.5) = - (-0.5 - 0.5) = 1 

H({B, C, D})= - (0.18 * log 2 0.18 + 0.82 * log 2 0.82) = -(-0.44 - 0.07) = 0.51 

5 Set {A, B, C} is thus seen to provide a more uniform coverage of volume 1600 than 
set {B, C, D}. 

In an additional preferred embodiment, the greedy algorithm is used to 
incrementally construct larger and larger sets, until a set containing the required 
number of views is obtained. As described for Figs. 10-11 above, the algorithm starts 

10 with an initial set of views (which may be the empty set), and in each iteration adds 
the view that yields the maximum improvement of the set score, until the required 
optimal set of views is obtained. 

In a further preferred embodiment of a method for the selection of an optimal 
set, multiple view sets are first formed from one or more scoring functions, and then a 

15 final selection is made of one of the resulting sets, substantially as described for Fig. 
12 above. Preferably the viewing parameters consist of at least one the following: 
detector location, detector orientation, viewing angle, material, thickness, collimator 
length, septa thickness, cell size, detection duration, time of detection, and a type of 
radiopharmaceutical. Knowledge about the conditions under which a particular set of 

20 measurements will be taken may be used to determine the optimal final set of views 
for a given data collection scenario. 

The abovedescribed methods may each be embodied as a computer program 
stored on a computer-readable storage medium. In the preferred embodiment, 
computer-readable storage medium contains a set of instructions for defining views 

25 for radioactive-emission measurements of the body structure. An analysis routine 
analyzes the radioactive-emission measurements obtained from a radioactive- 
emission-measuring probe, and a view definition routine dynamically defines further 
views for measurements, based on the analyzing. 

The active vision embodiments discussed above may be provided as 

30 modifications to improve reconstruction reliability. Reference is now made to Fig. 

17, which is a simplified flowchart of a method for stabilizing the reconstruction of an 
imaged volume with dynamic view definition, according to a preferred embodiment 
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of the present invention. Fig. 17 presents an embodiment in which the modifications 
to improve reconstruction reliability are performed by the active vision techniques 
described herein. Instep 1700 radioactive emission measurements are performed. In 
step 1710, the reliability of reconstruction is analyzed, for example utilizing the SVD 
analysis described above. In step 1720, a determination is made whether 
measurements are finished, preferably based on the reliability analysis performed in 
step 1710. If it is determined that modifications are required, further views are 
defined without ceasing the measurement process, and preferably utilizing active 
vision techniques. Further measurements are then performed at the newly-defined 
views, until an adequate reconstruction reliability is achieved or it is decided to 
tenninate the measurement process for other reasons. 

The above-described improved reconstruction reliability and active vision 
techniques enable high-quality reconstruction based on data collected from a limited 
collection of views. The ability to improve reconstruction reliability, both before and 
during imaging, reduces the need for repeated imaging to obtain additional 
measurements. Reconstructing the intensity distribution from a smaller quantity of 
collected data utilizing a stable reconstruction process simplifies the computational 
process. Furthermore, these methods enable resolving the current conflict between 
the relatively large-pixel detectors needed for measurement speed and data processing 
considerations, with the small-pixel detectors needed until now to obtain a high- 
resolution reconstruction. The abovedescribed embodiments are particularly suitable 
for medical imaging purposes, where a high-resolution image is needed and it is 
desired to minimize the difficulties of the patient undergoing the diagnostic testing or 
treatment. 
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• US Application 60/750,287, filed December 13, 2005: 

• US Application 60/750,334, filed December 15, 2005; 

• US Application 60/750,597, filed December 1 5, 2005; 

• US Application 60/800,845, filed May 17, 2006; 

• US Application 60/800,846, filed May 17, 2006; 

• Israel Application 171346, filed October 10, 2005; 

• Israel Application 1 72349, filed November 27, 2005 ; 

• US Application 60/741 ,440, filed December 2, 2005; 

. International Application PCT/IL2006/000059, filed January 15, 2006; 

• US Application 60/763,458, filed January 31, 2006; 

• International Application PCT/IL2006/000562, filed May 1 1 , 2006; 

• US Application 60/799,688, filed May 1 1 , 2006; and 

• US Application 60/8 1 6,970, filed June 28, 2006. 
Information of all of which is herein incorporated in entirety by reference. 

This application further incorporates by reference all the information of the 
International Application entitled "Imaging Protocols" which is being co-filed by the 
same assignee of the present invention on July 19, 2006. 

It is expected that during the life of this patent many relevant detection probes, 
detector types, radiation-based detection systems, algorithms for reconstruction and 
analyzing the reliability of reconstruction, and algorithms will be developed and the 
scope of the corresponding terms are intended to include all such new technologies a 
priori. 

It is appreciated that certain features of the invention, which are, for clarity, 
described in the context of separate embodiments, may also be provided in 
combination in a single embodiment. Conversely, various features of the invention, 
which are, for brevity, described in the context of a single embodiment, may also be 
provided separately or in any suitable subcombination. 
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Although the invention has been described in conjunction with specific 
embodiments thereof, it is evident that many alternatives, modifications and variations 
will be apparent to those skilled in the art. Accordingly, it is intended to embrace all 
such alternatives, modifications and variations that fall within the spirit and broad 
5 scope of the appended claims. All publications, patents and patent applications 
mentioned in this specification are herein incorporated in their entirety by reference 
into the specification, to the same extent as if each individual publication, patent or 
patent application was specifically and individually indicated to be incorporated 
herein by reference. In addition, citation or identification of any reference in this 
10 application shall not be construed as an admission that such reference is available as 
prior art to the present invention. 
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WHAT IS CLAIMED IS: 

1 . A method for stabilizing the reconstruction of an imaged volume, 
comprising: 

performing an analysis of the reliability of reconstruction of a radioactive- 
emission density distribution of said volume from radiation detected over a specified 
set of views; and 

defining modifications to at least one of a reconstruction process and a data 
collection process to improve said reliability of reconstruction, in accordance with 
said analysis. 

2. A method according to claim 1, further comprising calculating a measure of 
said reliability of reconstruction, said measure of reliability of reconstruction being 
for determining a necessity of performing said modifications. 

3. A method according to claim 1, further comprising: 

providing a detection probability matrix defining a respective detection 
probability distribution of said volume for each of said views; calculating the singular 
values of said detection probability matrix; 

identifying singular values as destabilizing singular values. 

4. A method according to claim 3, further comprising performing singular 
value decomposition (SVD) of said probability matrix. 

5. A method according to claim 3, wherein said identifying comprises: 

for each of said singular values, calculating a respective ratio of said singular 
value to a largest singular value; and 

specifying a singular value having a respective ratio below a specified 

threshold as destabilizing. 

6. A method according to claim 3, further comprising calculating a condition 
number of said probability matrix as a measure of said reliability of reconstruction. 
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7. A method according to claim 6, further comprising performing said 
identifying if said condition number exceeds a specified magnitude. 

8. A method according to claim 3, wherein said modification comprises 
defining constraints on said reconstruction process to reduce an effect of said 
destabilizing singular values upon a reconstruction of said volume. 

9. A method according to claim 3, further comprising identifying, for a 
destabilizing singular value, an associated unreliable voxel of said reconstructed 
volume. 

10. A method according to claim 3, further comprising identifying, for a 
destabilizing singular value, an associated unreliable linear combination of voxels of a 
reconstruction of said volume. 

1 1. A method according to claim 1, wherein said modification comprises 
smoothing a selected portion of said volume during said reconstruction. 

12. A method according to claim 1, wherein said modification comprises 
uniting multiple voxels of said volume into a single element of said reconstruction. 

13. A method according to claim 1, wherein said modification comprises 
defining views for imaging said volume. 

14. A method according to claim 13, wherein said defining is performed to 
obtain a required resolution. 

15. A method according to claim 13, wherein said defining comprises 
specifying multiple scan densities for imaging said volume. 
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16. A method according to claim 15, wherein each of said scan densities is for 
a different respective portion of said volume. 

17. A method according to claim 13, wherein said defining is performed 
dynamically during radioactive-emission measurements of said volume. 

1 8. A method according to claim 13, further comprising performing 
radioactive-emission measurements of the body structure, at a predetermined set of 
views; 

and wherein said modification comprises dynamically-defined views for 
continued imaging of said volume. 

19. A method according to claim 9, wherein said modification comprises 
defining a non-uniform scan for data collection. 

20. The method of claim 1, further comprising providing at least one intensity 
distribution, wherein an intensity distribution comprises a representation of a 
radioactive-emission density distribution of said volume. 

21. The method of claim 20, further comprising determining said intensity 
distribution from said reconstruction process. 

22. The method of claim 20, further comprising providing an omittance model 
as said intensity distribution. 

23. The method of claim 22, wherein said providing at least one emittance 
model includes modeling a radioactive-emission density distribution to form at least 
one modeled organ target. 

24. The method of claim 20, further comprising calculating an information- 
theoretic Fisher information as a measure of said reliability of reconstruction. 
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25. The method of claim 1, further comprising iteratively reconstructing a 
radioactive-emission density distribution of said volume. 

26. The method of claim 25, wherein said modification comprises performing 
said reconstructing using truncated SVD. 

27. The method of claim 25, wherein said reconstructing is by EM estimation. 

28. A reconstruction stabilizer, for improving the reliability of reconstruction 
of an imaged volume, comprising: 

a reliability analyzer, configured for performing an analysis of the reliability 
of reconstruction of a radioactive-emission density distribution of said volume from 
radiation detected over a specified set of views; and 

a modifier associated with said reliability analyzer, configured for defining 
modifications to at least one of a reconstruction process and a data collection process 
to improve said reliability of reconstruction, in accordance with said analysis. 

29. A reconstruction stabilizer according to claim 28, wherein said modifier is 
operable to define at least one view for imaging a volume. 

30. A reconstruction stabilizer according to claim 29, wherein said modifier is 
operable to define said at least one further view dynamically during radioactive- 
emission measurements of said volume. 

3 1. A reconstruction stabilizer according to claim 28, wherein said modifier is 
operable to apply at least one constraint to the reconstruction process 

32. A reconstruction stabilizer according to claim 31, wherein said analysis of 
the reliability of reconstruction comprises defining singular values of said detection 
probability matrix as destabilizing, and wherein said modifier is operable to apply a 
constraint to the reconstruction process to reduce an effect of said destabilizing 
singular values. 
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33. A reconstruction stabilizer according to claim 28, wherein said modifier is 
operable to designate a portion of said volume for smoothing during said 
reconstruction. 

34. A reconstruction stabilizer according to claim 28, wherein said modifier is 
operable to designate multiple voxels of said volume for uniting into a single element 
of said reconstruction. 

35. A system for generating a three-dimensional image of volume, the 
system comprising: 

a radiological imaging camera comprising a plurality of detectors configured 
for independent movement during data acquisition, and configured for detecting 
radiation emitted from said volume thereby to provide radiation data; 

a reconstructor, configured for performing an analysis of said radiation data so 
as to reconstruct a three-dimensional image of said volume; and 

a reconstruction stabilizer associated with said camera and said reconstructor, 
comprising: 

a reliability analyzer, configured for performing an analysis of the 
reliability of reconstruction of a radioactive-emission density 
distribution of said volume from radiation detected over a specified set 
of views; and 

a modifier, for defining modifications to at least one of a 
reconstruction process and a data collection process to improve said 
reliability of reconstruction, in accordance with said analysis, and for 
providing said modifications to at least one of said camera and said 
reconstructor. 

36. A method of radioactive-emission measurements of a body structure, 
comprising: 

performing radioactive-emission measurements of the body structure, at a 
predetermined set of views; 
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analyzing the radioactive-emission measurements; and 

dynamically defining further views for measurements, based on the analyzing. 

37. The method of claim 36, wherein said performing is in vivo. 

38. The method of claim 36, wherein said performing is for diagnostic 
purposes. 

39. The method of claim 36, wherein each of said views is associated with 
viewing parameters selected from the group consisting of: detector unit location, 
detector unit orientation, collection angle, and measurement duration. 

40. A method according to claim 36, wherein said analyzing comprises 
determining that a photon count at a given view yields a measurement error below a 
specified error value. 

4 1 . A method according to claim 40, wherein said defining comprises 
extending a duration of a current view to obtain a required error value. 

42. A method according to claim 36, wherein said analyzing comprises 
determining that a photon collection rate at a given view is within a specified range. 

43. A method according to claim 42, wherein said defining comprises 
selecting views to avoid regions of low interest. 

44. A method according to claim 42, wherein said defining comprises 
selecting views to concentrate on regions of high interest. 

45. A method according to claim 36, wherein said analyzing comprises 
identifying detector saturation at a given view. 
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46. A method according to claim 45, wherein said defining comprises 
selecting views to reinforce views that have saturated. 

47. A method according to claim 36, wherein said analyzing comprises 
forming a reconstruction of a radioactive-emission density distribution of said body 
structure. 

48. A method according to claim 47, wherein said analyzing comprises 
identifying textural edges within said reconstruction. 

49. A method according to claim 48, wherein said defining comprises 
selecting views at an angle to said textural edges. 

50. A method according to claim 49, wherein said angle is a substantially 
shaip angle. 

51. A method according to claim 47, wherein said analyzing comprises 
identifying volumetric boundaries within said reconstruction. 

52. A method according to claim 51, wherein said defining comprises 
selecting views at an angle to said volumetric boundaries. 

53. A method according to claim 52, wherein said angle is a substantially 
sharp angle. 

54. A method according to claim 47, wherein said defining comprises 
selecting views in close proximity to a suspected organ target. 

55. The method of claim 36, wherein said method is performed iteratively. 



56. The method of claim 47, wherein said method is performed iteratively. 
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57. A method according to claim 56, wherein said analyzing comprises 
determining a variance of a plurality of reconstructions. 

58. A method according to claim 57, wherein said defining comprises 
selecting views to reduce the variance of a region of high interest to a required level. 

59. A method according to claim 47, wherein said analyzing comprises 
determining an accuracy of said reconstruction. 

60. A method according to claim 47, wherein said analyzing comprises 
determining a resolution of said reconstruction. 

61. A method according to claim 36, wherein said dynamically defining 
comprises: 

providing a quality function for rating a quality of information by said quality 
function; and 

selecting a set of views to maximize said quality function. 

62. A method according to claim 61, wherein the quality function is based on 
an information theoretic measure of separability. 

63. A method according to claim 61, wherein the quality function is based on 
an information-theoretic Fisher information measure, for ensuring reliable 
reconstruction of a radioactive-emission density distribution. 

64. A method according to claim 61, further comprising increasing said 
selected set of views by at least one randomly selected view. 

65. The method of claim 61, wherein said dynamically defining comprises: 
providing a collection of views of said body structure; 

establishing said predetermined set of views as a current set of views; and 
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iteratively expanding said current set until a predefined number of further 
views are obtained, wherein a set expansion iteration comprises: 

forming a respective expanded set for each view in said collection of 
views and not in said current set, wherein an expanded set comprises 
said current set and a respective view; 

calculating a respective score for each of said expanded sets using said 
quality function; 

selecting said respective view of a highest-scoring expanded set as a 
further view; and 

equating said current set to said highest-scoring expanded set. 

66. The method of claim 36, further comprising selecting an initial 
predetermined set of views, selected by a method of: 

providing a model of a volume; 

providing a collection of views of said volume; 

providing a scoring function, by which any set of at least one view from said 
collection is scorable with a score that rates information obtained from said volume 
by said set; 

forming sets of views and scoring them, by said scoring function; and 
selecting said initial predetermined set of views from said collection, based on 
said scoring function. 

67. The method of claim 66, wherein each of said views is associated with 
viewing parameters comprising a detector location and orientation. 

68. The method of claim 66, wherein each of said views is associated with 
viewing parameters selected from the group consisting of: detector material, detector 
thickness, septa thickness, and collection angle. 

69. The method of claim 66, wherein each of said views is associated with 
viewing parameters selected from the group consisting: of measurement duration, 
time elapsed from the administration of the pharmaceutical to the measurement, 
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radiopharmaceutical half life, radioactive emission type, and radioactive emission 
energy. 

70. The method of claim 69, further comprising calculating for each of said 
views a respective detection probability distribution of said volume for each of said 
views, in accordance with respective viewing parameters and a volume attenuation 
coefficient. 

71. The method of claim 70, wherein said volume attenuation coefficient 
varies over said volume. 

72. The method of claim 66, wherein said selecting is to obtain a minimal-size 
set which attains a predefined score. 

73. The method of claim 66, wherein said selecting is to obtain a highest 
scoring set having a predefined number of views. 

74. The method of claim 66, wherein said collection of views comprises a 
quantized continuum of views. 

75. The method of claim 66, wherein said scoring function comprises an 
information-theoretic entropy measure, for ensuring uniform coverage of said volume. 

76. The method of claim 66, wherein said scoring function comprises a worst- 
case effectiveness for said given view over said volume. 

77. The method of claim 66, wherein said scoring function comprises an 
average effectiveness for said given view over said volume. 

78. The method of claim 66, further comprising providing at least one 
emittance model, wherein an emittance model comprises a representation of a 
radioactive-emission density distribution of said volume. 
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79. The method of claim 78, wherein said providing at least one emittance 
model includes modeling a radioactive-emission density distribution to form at least 
one modeled organ target. 

80. The method of claim 78, wherein said scoring function comprises an 
information-theoretic Fisher information measure, for ensuring reliable reconstruction 
of a radioactive-emission density distribution. 

8 1 . The method of claim 78, wherein: 

said providing at least one emittance model comprises providing a plurality of 
emittance models of substantially identical volumes, but different modeled organ 
targets; 

the forming sets of views from the collection of views and scoring them, 
comprises forming substantially identical sets of views for all the models and scoring 
the sets with respect to each model, based on the information theoretic measure of 
reliability; and 

the selecting one of the sets of views, based on its score, comprises selecting 
based on the average score for the plurality of models. 

82. The method of claim 78, wherein said scoring function comprises an 
information-theoretic likelihood measure, for ensuring separable reconstructions of 
said emittance models. 

83. The method of claim 78, wherein: 

said providing at least one emittance model comprises providing a pair of 
models of substantially identical volumes, but different modeled organ targets, 
wherein the difference between the modeled organ targets is defined by at least one 
delta; 

the scoring function is based on an information theoretic measure of 
separability; 
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the forming the sets of views from the collection of views and scoring them, 
comprises forming substantially identical sets of views for the pair of models and 
scoring the sets with respect to the pair; and 

the selecting one of the sets of views, based on its score, comprises selecting 

based on the score for the pair. 

84. The method of claim 78, wherein: 

said providing at least one emittance model comprises providing a plurality of 
pairs of models of substantially identical volumes, but different modeled organ 
targets, wherein the difference between the modeled organ targets for each pair is 

defined by at least one delta; 

the scoring function is based on an information theoretic measure of 

separability; 

the forming the sets of views from the collection of views and scoring them, 
comprises forming substantially identical sets of views for the plurality of pairs of 
models and scoring the sets with respect to each of the pairs; and 

the selecting one of the sets of views, based on its score, comprises selecting 
based on an average score for the plurality of pairs. 

85. The method of claim 78, wherein: 

the modeling the body-structure comprises providing at least one model of no 
radiation emission density distribution and at least two models of different radiation 
emission density distributions, so as to form different modeled organ targets for those 
two models; and 

the providing the scoring function comprises providing a scoring function as a 
combination of uniformity, separability and reliability. 

86. The method of claim 78, wherein said scoring function comprises a worst- 
case effectiveness for said given view over said set of emittance models. 



87. The method of claim 78, wherein said scoring function comprises an 
average effectiveness for said given view over said set of emittance models. 
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88. The method of claim 78, wherein said set of emittance models comprises 
a single emittance model. 

89. The method of claim 78, wherein said at least one emittance model 
corresponds to said body structure. 

90. The method of claim 89, wherein said body structure comprises one of a 
group of body structures comprising: a prostate, a heart, a brain, a breast, a uterus, an 
ovary, a liver, a kidney, a stomach, a colon, a small intestine, an oral cavity, a throat, a 
gland, a lymph node, the skin, another body organ, a limb, a bone, another part of the 
body, and a whole body. 

91. The method of claim 66, wherein said selecting comprises: 

initially establishing a current set of views from said collection of views; and 
iteratively expanding said current set until a predefined number of views are 
obtained, wherein a set expansion iteration comprises: 

forming a respective expanded set for each view not in said current set, 
wherein an expanded set comprises said current set and a respective 
view; 

calculating a respective score for each of said expanded sets using said 
scoring function; and 

equating said current set equal to a highest-scoring expanded set. 

92. The method of claim 91, wherein said initially established current set is an 
empty set. 

93. The method of claim 36, further comprising performing radioactive- 
emission measurements of the body structure, at said further views. 

94. A measurement unit for performing radioactive-emission measurements 
of a body structure, comprising: 
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a probe for performing the radioactive-emission measurements of the body 
structure, said probe being controllable to perform said measurements at a 
predetermined set of views; 

an analysis unit for analyzing the radioactive-emission measurements; and 
a view definer for dynamically defining further views for measurements, 
based on the analyzing. 

95 . A computer-readable storage medium containing a set of instructions for 
defining views for radioactive-emission measurements of the body structure, the set of 
instructions comprising: 

an analysis routine for analyzing the radioactive-emission measurements; and 
a view definition routine for dynamically defining further views for 
measurements, based on the analyzing. 
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